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INTRODUCTION 


Recently  the  AFFDL/Nielsen  subsonic  store  separation 
program.  Reference  1,  has  been  used  to  make  predictions  for 
comparisons  with  flight  and  wind-tunnel  captive  store  loads 
on  a  MK-83  bomb  on  an  F-4C  aircraft.  The  bomb  was  mounted 
on  the  bottom  station  of  a  triple  ejection  rack  (TER)  with 
dummy  bombs  mounted  on  the  two  shoulder  stations.  The  com¬ 
parisons  are  presented  in  Reference  2  and  indicate  that 
deficiencies  may  exist  in  the  TER  model  in  the  computer 
program . 

The  work  documented  in  this  report  had  two  main  objec¬ 
tives.  The  first  was  to  provide  a  data  base  which  could  be 
used  to  determine  where  deficiencies  in  the  computer  program 
exist.  The  second  objective  was  to  attempt  to  identify  the 
deficiencies  by  making  comparisons  of  computer  program 
results  with  the  data. 

The  test  program  was  conducted  in  the  4-Foot  Transonic 
Wind  Tunnel  ( 4T )  of  the  Propulsion  Wind  Tunnel  Facility 
(PWT)  at  the  Arnold  Engineering  Development  Center  (AEDC). 
The  testing  period  was  November  12  through  November  21, 

1979.  Reference  3  is  the  AEDC  document  describing  the  test 
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proqram.  The  test  program  was  jointly  sponsored  by  the 
Naval  Weapons  Center  (NWC)  and  the  Air  Force  Flight  Dynamics 
Laboratory  (AFFDL/FGC) . 

During  the  test  program,  force  and  moment  data  were 
obtained  using  both  captive  trajectory  system  (CTS)  and 
bracket  supported  store  models.  Flow-field  data  were 
obtained  using  a  20°  half-angle  conical  probe.  The 
testing  was  performed  at  Mach  numbers  from  0.6  to  0.95  and 
at  angles  of  attack  up  to  16°.  For  all  types  of  testing 
the  F-4C  parent  aircraft  configuration  was  built  up 
component  by  component  in  order  to  isolate  interference 
effects.  The  data  obtained  during  the  program  have  been 
recorded  on  magnetic  tape  and  FORTRAN  programs  have  been 
written  for  use  in  retrieving  the  data.  The  tapes  and 
programs  are  available  upon  request  from  Mr.  R.  E.  Smith, 
Code  3243,  Naval  Weapons  Center,  China  Lake,  CA  93555. 

The  next  section  of  this  report  will  expand  further  on 
the  purpose  and  scope  of  the  test  program,  describe  the  test 
apparatus,  and  present  the  test  conditions.  This  will  be 
followed  by  a  series  of  sub-sections,  each  one  devoted  to 
describing  a  specific  type  of  data.  Information  is 
presented  which  will  alLow  the  user  to  determine  the  AEDC 
"run"  number  for  the  data  in  which  he  is  interested.  These 
numbers  are  used  in  conjunction  with  the  magnetic  tapes  and 
computer  programs  to  retrieve  the  data.  Instructions  for 
doing  this  are  presented  and  the  computer  output  is 
described  . 

The  section  of  this  report  following  the  discussion  of 
the  experimental  program  presents  a  discussion  of  the 
experimental  results.  The  compatibility  of  the  captive  or 
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attached  loads  with  the  grid  loads  is  discussed.  Then 
phenomena  exhibited  by  the  data  during  the  parent  aircraft 
build-up  sequence  are  discussed  under  the  headings 

(a)  Attached  loads 

(b)  Grid  loads 

(c)  Flow  fields 

The  next  section  of  the  report  presents  comparisons 
between  the  <  xperimental  data  and  predictions  made  using  the 
computer  pr^jram  of  Reference  1.  This  section  will  try  to 
determine  limitations  of  the  computer  program  by  comparisons 
with  the  significant  effects  found  experimentally  as 
discussed  in  the  previous  section.  The  data  and  theory 
comparisons  will  be  used  to  determine  ways  of  improving  the 
program . 

The  last  section  of  the  report  will  present  specific 
recommendations  for  computer  program  modifications  which 
should  improve  the  accuracy  of  the  program. 


WIND-TUNNEL  TEST  PROGRAM 


PURPOSE  AND  SCOPE  OF  THE  TEST  PROGRAM 


The  purpose  of  the  test  program  was  to  provide  a  system¬ 
atic  set  of  data  which  could  be  used  in  evaluating  and 
improving  the  capability  of  the  computer  program  of 
Reference  1  to  predict  the  aerodynamic  forces  and  moments 
acting  on  stores  carried  in  a  TER  configuration.  The 
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capability  of  calculating  store  forces  and  moments  with  the 
store  in  the  carriage  position  is  of  interest  as  well  as 
the  capability  in  the  detached  position.  For  this  reason 
force  and  moment  tests  were  performed  using  both  bracket 
supported  and  CTS  supported  force  models.  In  order  to 
calculate  analytically  the  forces  and  moments  acting  on  a 
store,  the  aircraft  induced  flow  field  in  which  the  store 
is  immersed  must  be  predicted.  Flow-field  survey  data  were 
taken  to  provide  information  on  the  aircraft  induced  flow 
field . 

The  scope  of  the  program  was  to  obtain  these  three  types 
of  data  for  a  range  of  Mach  numbers,  0.6  to  0.95;  angles  of 
attack,  -4°  to  16°;  and  parent  aircraft  configurations.  The 
parent  aircraft  used  was  the  F-4C .  It  was  selected  because 
of  the  work  documented  in  Reference  2  which  initiated  this 
investigation.  The  parent  aircraft  was  built  up  component 
by  component  from  a  clean  F-4C  to  an  F-4C  with  wing  pylons, 
fuselage  pylon,  TER  on  the  inboard  wing  pylon,  and  stores  on 
the  shoulder  locations  of  the  TER.  The  store  used  in  the 
present  tests  was  the  Mk-83  bomb  used  in  the  work  of 
Reference  2.  The  component  by  component  build  up  was 
performed  so  that  the  effect  of  each  component  could  be 
determined.  For  example,  the  F-4C  with  inboard  left  wing 
pylon  and  TER  was  tested  with  and  without  Mk-83  stores  on 
the  shoulder  locations.  Differences  in  the  data  show  the 
effects  of  the  shoulder  stores  on  the  flow  field  and  the 
bottom  store  forces  and  moments. 
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TEST  APPARATUS 


Wind  Tunnel 


The  Aerodynamic  Wind  Tunnel  (4T)  is  a  closed-loop, 
continuous  flow,  variable-density  tunnel  in  which  the  Mach 
number  can  be  varied  from  0.1  to  1.3  and  can  be  set  at 
discrete  Mach  numbers  of  1.6  and  2.0  by  placing  nozzle 
inserts  over  the  permanent  sonic  nozzle.  The  nominal  range 
of  the  stagnation  pressure  is  from  400  to  3,400  pounds  per 
square  foot,  absolute.  The  test  section  is  4  feet  square 
and  12.5  feet  long  with  perforated,  variable  porosity  (0.5 
to  10  percent  open)  walls.  It  is  completely  enclosed  in  a 
plenum  chamber  from  which  the  air  can  be  evacuated,  allowing 
part  of  the  tunnel  airflow  to  be  removed  through  the 
perforated  walls  of  the  test  section.  A  more  complete 
description  of  the  test  facility  may  be  found  in  Reference  4. 


Model  Suppo rt  Systems 

For  this  test,  two  separate  and  independent  support 
systems  were  used.  The  aircraft  model  was  installed  inverted 
in  the  test  section  and  was  supported  by  an  offset  sting 
attached  to  the  main  pitch  sector.  For  loads  test¬ 

ing,  the  store  model  was  mounted  on  it  balance  fastened  to 
the  bottom  station  of  the  TER  on  the  aircraft  model.  For 
grid  aerodynamic  loads  and  flow-field  testing  the  store  model 
or  the  flow-field  probe  was  mounted  on  the  captive  trajectory 
support  ( CTS ) .  The  aircraft  model  was  removed  when  obtaining 
free-stream  data,  and  the  CTS  was  moved  upward  and  down¬ 
stream  in  the  tunnel  during  the  captive  loads  phase  to 
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minimize  interference.  Isometric  drawings  of  typical 
captive  loads,  grid  survey,  and  flow-field  testing  instal¬ 
lations  are  shown  in  Figure  1  of  Reference  3  along  with 
block  diagrams  of  the  computer  control  loops.  A  schematic 
showing  the  test  section  details  and  the  location  of  the 
models  in  the  tunnel  when  using  the  CTS  is  shown  in  Figure  2 
of  that  reference.  A  further  description  of  the  CTS  rig  can 
be  found  in  Reference  4. 


Parent  Aircraft  Model  Components  and  Stores 

The  basic  details  of  the  five-percent  (l/20th)  scale 
model  of  the  F-4C  are  shown  in  Figure  1.  The  model  is 
geometrically  similar  to  the  full-scale  aircraft  except  that 
the  part  aft  of  the  engine  exhausts  has  been  removed  to 
minimize  CTS  interference.  This  removal  does  not  influence 
the  results  of  this  test  progi-am.  The  F-4C  model  has  flow¬ 
through  engine  inlets  with  subsonic  and  supersonic  exhaust 
chokes.  The  subsonic  chokes  were  used  during  the  test  pro¬ 
gram.  During  part  of  the  captive  loads  phase  of  the  test 
program  a  total  pressure  rake,  containing  13  total  pressure 
orifices,  was  mounted  just  aft  of  the  right-side  engine 
exhaust  choke.  The  interior  surface  of  the  choke  itself  was 
instrumented  with  six  static  pressure  orifices.  These 
pressure  measurements  were  made  in  order  that  a 
calculation  of  the  ratio  of  the  inlet  velocity  relative  to 
the  free-stream  velocity  could  be  made.  This  velocity  ratio 
influences  the  input  data  to  the  computer  program  of 
Reference  1.  Details  and  dimensions  of  the  total  pressure 
rake  and  the  static  pressure  instrumentation  of  the 
exhaust  choke  are  shown  in  Figure  6  of  Reference  3. 
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Details  of  the  F-4C  pylons  are  shown  in  Figure  2.  For 
most  tests  only  the  left  wing  inboard  pylon  was  used.  All 
testing  was  done  under  the  left  wing.  Pylons  were  not 
installed  on  the  right  wing  for  any  of  the  tests. 

The  triple  ejector  rack  (TER)  used  in  the  test  program 
is  shown  in  Figure  3.  The  TER  was  used  in  conjunction  with 
the  left-wing  inboard  pylon.  This  rack  model  simulates  sway- 
braces  and  ventilating  passages. 

The  stores  which  were  used  are  shown  in  Figure  4.  They 
are  models  of  the  Mk-83  bomb.  Figure  4(a)  shows  the  actual 
bomb  shape  and  Figure  4(b)  shows  a  model  which  was  modified 
for  CTS  sting  support.  Models  with  identical  body  shapes  but 
without  tail  fins  were  also  used.  The  dummy  stores  used  on 
the  shoulder  locations  of  the  TER  were  the  shape  shown  in 
Figure  4(a).  Both  finned  and  unfinned  models  were  used. 

Four  force  models  were  used  during  the  captive  force  and 
moment  tests.  Tests  were  performed  with  both  configurations 
shown  in  Figure  4  with  and  without  tail  fins. 


TEST  CONDITIONS 

The  Mach  number  range  used  in  the  test  program  was  0.6 
to  0.95.  All  testing  was  performed  at  a  nominal  Reynolds 

C 

number  per  foot  of  3.5X10  .  The  store  angle  of  attack  was 
varied  from  -4°  to  16°.  Since  the  store  is  oriented  on  the 
TER  one  degree  nose  down  relative  to  the  aircraft  reference 
waterline  (aircraft  angle  of  attack  is  measured  relative  to 
this  line),  the  aircraft  angle  of  attack  varied  from  -3°  to 
17°.  Nominal  values  of  other  tunnel  parameters  are  listed 
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FS 

10.652 

(a)  Actual  configuration. 


ALL  DIMENSIONS  IN  INCHES 


(b)  Configuration  modified  for  sting  support. 


FIGURE  4.  Mk-83  Bomb  Models. 
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in  the  following  table  as  a  function  of  the  Mach  numbers 
used  during  the  tests. 


TABLE  1.  Nominal  Test  Conditions. 


M 

00 

Pt 

00 

psfa 

Soo 

psf 

T 

t 

00 

0  R 

*e_ 

ft"1 

0.6 

2186 

432 

549.6 

3.5xl06 

0.7 

1975 

489 

0.8 

1835 

539 

0.9 

1741 

584 

0.95 

1708 

604 

i 

PARENT  AIRCRAFT  CONFIGURATIONS 


Eight  parent  aircraft  configurations  were  used  during 
the  test  program.  They  are  listed  in  Table  2.  The  first 
column  lists  the  eight  numbers.  These  numbers  will  be  used 
in  the  next  sections  of  this  report  which  describe  the 
tests.  The  remaining  eight  columns  list  model  components 
which  could  be  attached  to  the  parent  F-4C .  These 
components  are 


P 

P 

P 


I 

O 


fuselage  centerline  pylon 
left-wing  inboard  pylon 
left-wing  mu-board  pylon 
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TABLE  2.  Parent  Aircraft  Configurations. 


CONFIG. 
NO.  * 

P£ 

PI 

po 

■ 

(S2>U 

(S3}  U 

( S2 )  F 

(S3>F 

1 

2 

X 

3 

X 

X 

4 

X 

X 

X 

5 

X 

X 

X 

X 

6 

X 

X 

X 

X 

7 

X 

X 

X 

X 

X 

X 

8 

X 

X 

X 

* 

Configuration  composed  of  basic  F-4C  model  with 
addition  of  components  indicated  in  line  across 
table . 


i 


NWC  TP  6210 


T  TER  attached  to  left-wing  inboard  pylon 

( 5 2 )  U  unfinned  dummy  store  attached  to  outboard 

shoulder  station  on  TER 

( 5 3 )  u  unfinned  dummy  store  attached  to  inboard 
shoulder  station  on  TER 

finned  dummy  store  attached  to  outboard 
shoulder  station  on  TER  with  fins  vertical 
and  horizontal 

finned  dummy  store  attached  to  inboard 
shoulder  station  on  TER  with  fins  vertical 
and  horizontal 

The  dummy  store  numbering  system  is  shown  in  Figure  5.  A 
check  in  a  box  opposite  a  configuration  number  in  Table  2 
indicates  this  model  component  was  attached  to  the  clean 
F-4C  . 

CAPTIVE  FORCE  AND  MOMENT  TESTS 

The  captive  force  and  moment,  phase  of  the  test  program 
consisted  of  taking  six-component  force  and  moment  data  on 
four  different  Mk-83  bomb  models.  The  two  store  shapes 
shown  in  Figure  4  were  used  with  and  without  tail  fins.  The 
test  installation  is  shown  in  Figure  6.  A  strain-gage 
balance  was  installed  internally  in  the  store.  It  was 
supported  by  a  bracket  protruding  from  the  TER.  The  balance 
leads  were  run  through  the  TER,  pylon,  and  F-4C .  All 
testing  was  done  with  the  force  model  mounted  on  the  bottom 
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station  of  the  TER.  The  finned  stores  were  always  attached 
to  the  balance  with  the  fins  rolled  45°  from  the 
vertical  and  horizontal . 

The  captive  aerodynamic  loads  data  were  obtained  at 
aircraft  angles  of  attack  from  -3°  to  17°,  store  angles  of 
attack  from  -4°  to  16°,  at  zero  degrees  sideslip  angle.  All 
angle  of  attack  polars  were  run  automaticaly  utilizing 
online  computer  facilities  to  calculate  the  control  commands 
to  set  the  aircraft  model  attitude. 

The  force  and  moment  data  for  the  store  were  reduced  to 
coefficient  form  using  the  standard  4T  data  reduction  pro¬ 
grams.  The  forces  and  moments  are  in  the  store  body  axis 
system  with  the  moments  taken  about  the  store  center  of 
gravity  location  shown  in  Figure  4.  Positive  senses  of  the 
forces  and  moments  and  the  store  body  axis  system  are  shown 
in  Figure  7.  During  the  present  tests  the  fins  were  rolled 
45°  from  the  orientation  shown  in  the  figure.  The  estimated 
aerodynamic  coefficient  uncertainties  are  listed  in  Table  3. 
They  are  taken  froin  Reference  3. 


TABLE  3.  Aerodynamic  Coefficient  Uncertainties. 


M 

OO 

CN 

CY 

CA 

Cm 

C 

n 

0.60 

±0.03 

±0.03 

±0.04 

±0.04 

±0.03 

±0.01 

0.70 

±0.02 

±0.02 

±0.04 

±0.04 

±0.02 

±0.01 

0.80 

±0.02 

±0.02 

±0.03 

±0.03 

±0.02 

±0.01 

0.90 

±0.02 

±0.02 

±0.03 

±0.03 

±0.02 

±0.01 

0.95 

+  0.02 

+  0.02 

+  0.03 

±0.03 

±0.02 

±0.01 
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FIGURE  7.  Axis  System  for  Store  Showing  Positive 
Sense  of  Axes,  Forces  and  Moments,  and  Velocities. 
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During  the  first  part  of  the  captive  loads  phase  the  13 
pressures  measured  by  the  total  pressure  rake  just  aft  of 
the  right-side  engine  exhaust  choke  and  the  six  static 
pressures  measured  inside  the  choke  were  recorded.  These 
pressures  were  used  to  calculate  the  engine  duct  exit  Mach 
number,  exit  mass  flow  rate,  capture  ratio  (exit  mass  flow 
rate  divided  by  theoretical  (reference)  inlet  mass  flow 
rate),  and  exit  velocity  divided  by  free-stream  velocity. 

The  following  procedure  was  used  to  calculate  these 
quantities.  The  average  static  pressure,  p,  was  calculated 
as  was  the  average  total  pressure,  pt .  The  ratio  p/pt  was 
determined.  If  this  ratio  was  greater  than  0.5283,  the  exit 
flow  was  subsonic  and  the  exit  Mach  number  was  calculated 
from 


M 


e 


(1) 


For  all  tests  performed  during  this  test  program  the  flow 
was  subsonic  although  provision  was  made  in  the  data 
reduction  program  to  handle  supersonic  exit  flows, 
p/pt  <  0.5283. 


With  the  exit  Mach  number,  Me,  determined,  the  exit 
mass  flow  is  calculated. 


me  * 


■We11 


2  “3 

+  0.210 

e 


0.8843 


1/2 


(2) 
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the  capture  ratio  is  calculated. 


c_  =  m  / 
R  e 


p.  M  A. (1  + 


0.2M2)"3 


0.8843 


1/2 


(3) 


as  is  the  ratio  of  exit  velocity  to  free-stream  velocity. 


V 


M 

€ 


1  +  0.2M“ 


1  +  0.2M 


1/2 


(4) 


In  Equation  (2)  the  cross  sectional  area  of  the  duct 

_  O 

exit,  Ae,  used  was  7.4662x10  °  square  feet,  model  scale.  The 
inlet  area,  A^ ,  appears  in  Equation  (3).  In  the  data  reduc- 
tion  a  value  of  1.705x10  square  feet,  model  scale,  was  used. 
This  value  was  provided  to  AEDC  by  McDonald  Douglas  and  is 
the  value  they  use.  The  cross-sectional  area  of  the  inlet 

_  O 

opening  on  the  wind-tunnel  model  was  l.lx10  square  feet, 
model  scale. 


Summary  of  Tests 

The  captive  loads  tests  are  summarized  in  Table  4. 

Column  1  lists  the  parent  aircraft  configuration  number  from 
Table  2.  Column  2  lists  the  store  force  and  moment  model 
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TABLE  4.  Captive  Force  and  Moment  Tests. 


CONFIG 

NO. 


See  footnote  at  end  of  table. 
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TABLE  4.  (Contd.) 


CONFIG. 

NO. 

FORCE 

MODEL 

M 

CO 

1 - 5 - 

s 

'INITIAL 

DEG. 

a 

s 

FINAL 

DEG. 

Aa 

s 

DEG. 

AEDC 

RUN  NO. 

6 

mmtmm 

0.6 

- 

4 

16 

2 

153 

6 

0.95 

i 

154 

5 

0.6 

164 

5 

0.95 

165 

3 

. j 

0.6 

160 

3 

i 

0.95 

1 

161 

+Repeat  with  total  pressure  rake  aft  of  right  engine  exhaust 
removed. 

Repeat  with  pressure  rake  removed  and  hole  in  TER  station 
no.  1  pylon  filled. 
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attached  to  the  parent  aircraft  at  TER  station  no.  1.  The 
store  model  designations  are 

SAF  Mk-83  actual  configuration  force  model  with 

tail  fins  (see  Figure  4(a)) 

SAU  same  as  S but  without  tail  fins 

SMF  Mk-83  modified  configuration  force  model  with 

tail  fins  (see  Figure  4(b)) 

SMU  same  as  SMp  but  without  tail  fins 

The  third  column  lists  the  free-stream  Mach  number,  M^,  and 
the  next  three  columns  list  the  initial  store  angle  of 
attack,  the  final  store  angle  of  attack,  and  the  increment 
in  angle  of  attack.  The  parent  aircraft  angle  of  attack, 
measured  relative  to  waterline,  WL,  zero  shown  in  Figure  1 
is  one  degree  greater  than  the  store  angle  of  attack. 

The  last  column  lists  the  AEDC  run  number  for  the  data 
obtained  for  that  particular  configuration  and  Mach  number. 
This  column  provides  input  data  for  the  data  retrieval 
computer  program. 

Some  explanation  of  the  two  footnotes  on  the  table  is 
required.  During  run  numbers  86  through  101,  the  total 
pressure  rake  just  aft  of  the  right-side  engine  exhaust 
was  present.  The  runs  marked  with  the  superscript  +  were 
made  with  the  rake  removed  to  determine  the  effect  of  the 
rake  on  the  store  forces  and  moments.  The  rake  was  found 
not  to  effect  the  forces  and  moments.  It  was  also  deter¬ 
mined  as  the  testing  progressed  that  for  runs  86  through 
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122  a  small  screw  hole  in  the  leading  edge  of  the  TER 
station  number  1  pylon  was  not  filled  with  putty  or  wax. 
It  was  filled  and  the  runs  with  the  superscript  *  were 
made  to  determine  the  effect.  It  was  found  that  the  dif¬ 
ferences  in  the  forces  and  moments  with  and  without  the 
hole  filled  were  within  the  data  uncertainties  listed  in 
Table  3. 


Use  of  Data  Retrieval  Computer  Program 

The  captive  loads  force  and  moment  data  are  all 
written  on  the  first  file  of  a  magnetic  tape.  The  second 
file  on  the  tape  contains  the  CTS  or  grid  force  and  moment 
data  to  be  described  in  the  next  section  of  this  report. 

A  computer  program  has  been  written  which  allows  the 
captive  loads  data  for  run  numbers  listed  in  Table  4  to 
be  retrieved  from  the  tape  and  tabulated.  To  retrieve 
these  data  the  force  and  moment  data  tape  must  be  positioned 
at  the  beginning  of  the  first  file. 

A  listing  of  the  computer  program  is  presented  in 
Appendix  A.  The  first  card  of  input  is  in  a  3A4  format 
and  contains  the  word  "CAPTIVE"  beginning  in  column  1. 

This  is  used  to  compare  with  a  header  on  the  tape  to 
assure  that  the  correct  tape  is  mounted  and  that  it  is 
positioned  at  the  correct  file.  The  second  input  card  in 
an  15  format  specifies  the  number  of  run  numbers  of  data 
to  be  retrieved.  The  specific  run  numbers  from  Table  4 
are  input  on  the  next  cards,  eight  to  a  card  in  8F10.0 
format.  The  run  numbers  can  be  in  any  order.  The  com¬ 
puter  program  sorts  them  into  numerical  sequence  before 
searching  the  tape. 
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Description  of  Tabulated  Data 

There  is  one  page  of  output  for  each  run  number  of  data. 
A  sample  is  shown  in  Figure  8.  The  first  line  of  output 
identifies  the  type  of  data,  in  this  case  captive  loads;  the 
AEDC  report  number  documenting  this  test  (Reference  3);  and 
the  ARO  project  number  and  test  number.  The  next  two  lines 
list  the  run  number,  parent  aircraft  configuration,  and 
store  model.  The  following  four  lines  of  output  list  the 
free-stream  test  conditions.  The  remainder  of  the  page 
tabulates  the  data  obtained.  When  the  last  four  columns  are 
equal  to  0.000,  the  total  pressure  rake  just  aft  of  the 
right-side  engine  exhaust  was  not  used. 

The  nomenclature  used  in  the  tabulated  output  is  con¬ 
tained  in  the  following  list. 

ALPHA  PAR  angle  of  attack  of  the  parent  aircraft,  deg. 

ALPHA  STORE  angle  of  attack  of  the  store,  deg. 

BETA  STORE  sideslip  angle  of  the  store,  positive 
nose  to  the  right,  deg. 

CAT  store  axial-force  coefficient  uncorrected  for 

base  pressure,  axial  fore c/(C  S) 

CLL  store  rolling-moment  coefficient,  rolling 

moment/ (Q  S  L) 

CLM  store  pitching-moment  coefficient,  pitching 

moment/ (Q  S  L) 
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CLN  store  yawing-moment  coefficient,  yawing 

moment/ (Q  S  L) 

CN  store  normal-force  coefficient,  normal 

force/ (Q  S) 

CR  capture  ratio  calculated  using  Equation  (3) 

CY  store  side-force  coefficient,  side 

force/ ( Q  S ) 

L  reference  length  used  in  computing  moment 

coef f icients ,  0.7  inch  model  scale 

MACH  NO  free-stream  Mach  number 

MDOTN  duct  exit  mass-flow  rate  calculated  using 

Equation  (2),  lL>r../sec 

MNE  duct  exit  Mach  number 

NCP  normal-force  eenter-of-pressure  location, 

CLM/CN;  store  diameters,  L,  from  the 
center  of  gravity,  positive  ahead 

P  free-stream  static  pressure, 

PT  free-stream  total  pressure,  psfa 

Q  free-stream  dynamic  pressure,  psf 

RE/FT  free-stream  Reynolds  number  per  foot,  ft-'*' 
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S  reference  area  used  in  computing  force  and 

moment  coefficients,  0.385  square  inches  model 
scale 

T  free-stream  static  temperature,  °R 

TT  free-stream  total  temperature,  °R 

V  free-stream  velocity,  ft/sec 

VR  ratio  of  engine  exit  velocity  to  free-stream 

velocity  calculated  using  Equation  (4) 

YCP  side-force  center-of-pressure  location, 

CLN/CY;  store  diameters,  L,  from  the 
center  of  gravity,  positive  ahead 

GRID  FORCE  AND  MOMENT  TESTS 

The  grid  force  and  moment  phase  of  the  test  program 
consisted  of  taking  six-component  force  and  moment  data  on 
two  different  Mk-83  bomb  models.  The  store  shape  shown  in 
Figure  4(b)  was  used  with  and  without  tail  fins.  A  sketch 
of  the  test  installation  is  shown  in  Figure  9.  The  paieuc 
aircraft  is  mounted  on  tne  main  pitch  sector  and  the  store 
is  mounted  on  an  internal  balance.  A  sting  protruding  from 
the  base  of  the  store  attaches  the  store  and  balance  tc  the 
CTS .  The  finned  store  was  always  attached  to  the  balance 
with  the  fins  rolled  45°  from  the  vertical  and  horizontal. 
This  is  the  zero  degree,  4>s  =  0°,  roil  '~>r  i  entat  ion  of 
the  store. 
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Grid  force  and  moment  data  were  obtained  at  aircraft 
angles  of  attack  from  1°  to  17°,  store  angles  of  attack  from 
0°  to  16°,  at  zero  degrees  sideslip  angle.  Testing  was 
accomplished  by  placing  the  aircraft  at  a  particular  angle 
of  attack.  The  store  was  then  placed  as  close  as  possible 
to  the  carriage  position  at  TER  station  number  1 .  Data  were 
taken  at  a  series  of  points  as  the  store  was  moved  away  from 
this  position  parallel  to  itself. 

Free-stream  data  were  obtained  on  each  of  the  two 
stores,  alone.  For  these  tests  the  parent  aircraft  model 
was  removed  and  data  were  taken  for  store  angles  of  attack 
from  -4°  to  16°.  Data  were  obtained  on  the  finned  store 
through  a  roll  angle  range  of  0°  to  90°  in  22.5°  increments. 

The  force  and  moment  data  for  the  store  were  reduced  to 
coefficient  form  using  the  standard  4T  data  reduction  pro¬ 
grams.  The  forces  and  moments  are  in  the  store  body  axis 
system  and  the  moments  are  taken  about  the  store  center  of 
gravity  location  shown  in  Figure  4.  Positive  senses  of  the 
forces  and  moments  are  shown  in  Figure  7.  As  the  store  was 
rolled  during  the  free-stream  tests,  the  store  X,Y,Z 
coordinate  system  rolled  with  the  store.  The  aerodynamic 
coefficient  uncertainties  are  listed  in  Table  3.  They  are 
taken  from  Reference  3. 


Summary  of  Tests 

The  store  free-stream  tests  are  summarized  in  Table  5. 
Column  1  lists  the  store  force  model .  The  store  model 
designations  are 


< 
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TABLE  5.  Store  Free-Stream  Force 
and  Moment  Tests. 
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SMF  Mk-83  modified  configuration  force  model 

with  tail  fins  (see  Figure  4(b)) 

SMU  same  as  SMF  but  without  tail  fins 

The  second  column  gives  the  free-stream  Mach  number,  M^, 
and  the  third  column  the  store  roll  angle,  $  .  Positive 
roll  is  clockwise  when  viewed  from  the  rear.  The  next 
three  columns  list  the  initial  store  angle  of  attack,  the 
final  store  angle  of  attack,  and  the  increment  in  angle 
of  attack.  The  last  column  presents  the  AEDC  run  number 
for  the  data  obtained  for  that  particular  Mach  number,  roll 
angle  combination.  This  column  provides  input  data  for  the 
data  retrieval  computer  program. 

Table  6  summarizes  the  grid  force  and  moment  tests  with  a 
parent  aircraft  present.  Column  1  lists  the  parent  aircraft 
conf iguration  number  from  Table  2.  Column  2  indicates  which 
of  the  two  force  models  was  mounted  on  the  balance.  The 
free-stream  Mach  number,  Mot ,  is  listed  in  column  3  and  the 
store  angle  of  attack,  ag,  in  column  4.  The  parent  aircraft 
angle  of  attack  measured  with  respect  to  waterline,  WL,  zero 
shown  in  Figure  1  is  one  degree  greater  than  the  store  angle 
of  attack.  The  fifth  column  indicates  which  Zg/D  schedule 
was  used  for  that  particular  run.  D  is  the  store  maximum 
diameter,  0.7  inch  model  scale.  The  origin  of  the  Xg ,  Yg, 

Zg  coordinate  system  is  shown  in  Figure  5  and  is  located  at 
the  tip  of  the  store  nose  when  the  store  is  in  the  carriage 
position  on  the  bottom  station  of  the  TER.  The  two  Zg/D 
schedules  are  shown  in  the  last  two  columns  of  the  first 
page  of  Table  6.  The  value  of  Zg/D  =  0.0  could  not  be 
obtained  since  this  would  result  in  the  store's  touching  the 
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TABLE  6.  Grid  Force  and  Moment  Tests. 


CONFIG. 

NO. 


FORCE 

MODEL 


M 


s 

DEG. 


Zs/D 

SCHD. 


AEDC 
RUN  NO. 


Z  g/D  SCHDS. 


ZA 


ZB 


1  s 

0, 

.  6 

0 

ZA 

343 

0.0 

0.0 

MF 

4 

344 

.  4 

.1 

8 

345 

.  8 

.  2 

12 

346 

1.2 

.  3 

16 

347 

1.6 

.4 

0 

.95 

0 

352 

2.0 

.  6 

4 

351 

2.4 

.  8 

8 

350 

2.8 

1.0 

12 

349 

3.2 

1.2 

16 

348 

3.6 

1.4 

2 

0. 

.  6 

0 

338 

4.0 

1.6 

J 

L 

y 

L 

4 

339 

4.4 

2.0 

1 

f 

8 

f 

340 

4.8 

2.4 

3 

0. 

.  6 

0 

ZB 

300 

5.2 

2.8 

4 

301 

5.6 

3.2 

8 

302 

6.0 

3.6 

12 

303 

4.0 

16 

304 

4.4 

0. 

.8 

0 

310 

4.8 

4 

309 

5.2 

8 

308 

5.6 

12 

307 

6.0 

16 

305 

0. 

.95 

0 

311 

4 

312 

8 

313 

12 

314 

16 

315 

< 

) 

0, 

.  6 

0 

249,327 

4 

250 

8 

251 

12 

254 

16 

255  . 

0. 

.  7 

0 

258,328 

J 

L 

4 

257 

8 

256  , 

0. 

8 

0 

260,329 

1 

4 

261 

1 

1 

r 

8 

\ 

262 

See 

footnote 

at 

end  of 

table . 
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CONFIG. 

NO. 


FORCE 

MODEL 


M 

00 

as 

VD 

DEG. 

SCHD. 

AEDC 
RUN  NO. 
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TABLE  6.  (Contd.) 


CONFIG. 

FORCE 

M 

a 

Z  /D 

AEDC 

NO. 

MODEL 

s 

DEG. 

SCHD. 

RUN  NO. 

ZA  ZB 

6  S„rT  0.95  0  ZB  245 


i 

*u  1  4 

244 

t 

f  8 

243 

3 

0.6  0 

232 

i  4 

233 

I  8 

234 

0.95  0 

237 

i  4  , 

236 

!  (  | 

(is' 

235 

(see  1st 
page  of 
table) 


* 

Repeat  with  TER  station  number  1  sway  braces 
removed . 
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parent  aircraft.  The  first  value  of  Zg/D  was  usually 
between  0.05  and  0.10. 

Column  6  of  Table  6  lists  the  AEDC  run  number  for  the 
data  obtained  for  that  particular  parent  aircraft  configura¬ 
tion,  Mach  number,  and  angle  of  attack.  This  column 
provides  input  data  for  the  data  retrieval  computer  program. 


Use  of  Data  Retrieval  Computer  Program 

The  grid  force  and  moment  data  are  all  written  on  the 
second  file  of  a  magnetic  tape.  The  first  file  on  the  tape 
contains  the  captive  loads  data  which  were  previously 
described.  A  computer  program  has  been  written  which  allows 
the  grid  data  for  run  numbers  listed  in  Tables  5  and  6  to  be 
retrieved  from  the  tape  and  tabulated.  To  retrieve  these 
data  the  force  and  moment  data  tape  must  be  positioned  at 
the  beginning  of  the  second  file. 

A  listing  of  the  computer  program  is  presented  in 
Appendix  B.  The  first  card  of  input  is  in  a  3A4  format  and 
contains  the  word  "GRID"  beginning  in  column  1.  This  is 
used  to  compare  with  a  header  on  the  tape  to  assure  that  the 
correct  tape  is  mounted  and  that  it  is  positioned  at  the 
correct  file.  The  second  input  card  in  an  15  format 
specifies  the  number  of  run  numbers  of  data  to  be  retrieved. 
The  specific  run  numbers  from  Tables  5  and  6  are  input  on 
the  next  cards,  eight  to  a  card  in  8F10.0  format.  The  run 
numbers  can  be  in  any  order.  The  computer  program  sorts 
them  into  numerical  sequence  before  searching  the  tape. 
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Description  of  Tabulated  Data 

There  is  one  page  of  output  for  each  run  number  of  data. 
A  sample  is  shown  in  Figure  10.  The  first  line  of  output 
identifies  the  type  of  data,  in  this  case  grid  force  data; 
the  AEDC  report  number  documenting  this  test  (Reference  3); 
and  the  ARO  project  number  and  test  number.  The  next  two 
lines  list  the  run  number,  parent  aircraft  configuration, 
and  store  model.  The  following  four  lines  of  output  list 
the  free-stream  test  conditions.  The  remainder  of  the  page 
tabulates  the  data  obtained. 

The  nomenclature  used  in  the  tabulated  output  is  con¬ 
tained  in  the  following  list. 

ALPHA  PAR  angle  of  attack  of  the  parent  aircraft,  deg. 

ALPHA  STORE  angle  of  attack  of  the  store,  deg. 

BETA  STORE  sideslip  angle  of  the  store,  positive  nose 
to  the  right,  deg. 

CAT  store  axial-force  coefficient  uncorrected  for 

base  pressure,  axial  force/(Q  S) 

CLL  store  rolling-moment  coefficient,  rolling 

moment/ (Q  S  L) 

CLM  store  pitching-moment  coefficient,  pitching 

moment/ (Q  S  L) 

Ci:  store  normal-force  coefficient,  normal 

force/(Q  S) 
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CY 

DPHI 


store  side-force  coefficient,  side  force/(Q  S) 

store  roll  angle,  positive  clockwise  looking 
upstream,  deg. 


DTHA  store  total  angle  of  attack,  angle  between 

store  longitudinal  axis  and  the  free-stream 
velocity  vector,  deg. 

L  reference  length  used  in  computing  moment 

coefficients,  0.7  inch  model  scale 


MACH  NO 


free-stream  Mach  number 


NCP  normal-force  center-of-pressure  location, 

CLM/CNr  store  diameters,  L,  from  the 
center  of  gravity,  positive  ahead 

P  free-stream  static  pressure,  psfa 


PT 


free-stream  total  pressure,  psfa 


Q 

RE/ FT 

S 


T 

TT 


free-stream  dynamic  pressure,  psf 

free-stream  Reynolds  number  per  foot,  ft-^ 

reference  area  used  in  computing  force  and 
moment  coefficients,  0.385  square  inches  model 
sea  1  e 

free-stream  static  temperature,  °R 
free-stream  total  temperature,  °R 
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free-stream  velocity,  ft/sec 


Xs-coordinate  of  the  store  nose  in  the  Xg ,  Yg , 
Z  coordinate  system  of  Figure  5,  in. 


side-force  center-of-pressure  location, 
CLN/CY;  store  diameters,  L,  from  the 
center  of  gravity,  positive  ahead 


Yg-coordinate  of  the  store  nose  in  the  Xg ,  Yg , 
Z  coordinate  system  of  Figure  5,  in. 


Zs-coordinate  of  the  store  nose  divided  by  the 
store  maximum  diameter  (0.7  in)  in  the  Xg ,  Yg , 
Zs  coordinate  system  of  Figure  5 


FLOW-FIELD  SURVEY  TESTS 

The  flow-field  survey  phase  of  the  test  program  con¬ 
sisted  of  measuring  the  flow  field  in  the  region  beneath  the 
parent  aircraft  which  the  bottom  store  on  the  TER  would 
occupy  in  and  near  the  carriage  position.  Details  of  the  20 
degree  half-angle  conically  tipped  probe  used  during  the 
tests  are  shown  in  Figure  11.  The  test  installation  is 
similar  to  that  shown  in  Figure  9  for  the  grid  force  and 
moment  tests  except  that  the  probe  is  mounted  on  the  CTS 
rather  than  the  balance  and  force  model.  The  p^obe  has  four 
equally  spaced  static  pressure  orifices  on  the  cone  surface 
and  a  total  pressure  orifice  at  the  cone  apex. 

Flow-field  data  were  taken  at  aircraft  angles  of  attack 
from  1°  to  17°.  The  probe  angle  of  attack  was  kept  the  same 
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as  the  bottom  store  on  the  TER,  that  is,  one  degree  less 
than  the  parent  aircraft.  The  origin  of  the  probe 
coordinate  system  is  the  Xp,  Yp,  Zp  system  shown  in  Figure 
5.  The  probe  longitudinal  axis  and  the  Xp  axis  are 
parallel . 

Data  were  obtained  by  making  axial  traverses  at  a  series 
of  Yp,  Zp  positions.  These  locations  are  shown  in  the 
following  sketch  by  the  crosses. 


TRAVERSE 

LOCATION 


0.35**— 


Y 


P 


t 


The  dashed  circle  represents  the  Mk-83  bomb  in  the  carriage 
position  on  the  TER.  The  Xp  range  for  most  of  the  traverses 
was  from  Xp  =  -0.15  inch  to  Xp  =  -5.85  inches,  approximately 
the  length  of  the  Mk-83  model.  For  two  of  the  traverses, 

Yp  =  0.0  and  Zp  =  0.0  and  0.70,  the  range  was  extended  to 
-0.85  inches. 
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Standard  4T  data  reduction  programs  were  used  to  calcu¬ 
late  the  local  velocity  components  and  flow  angles.  These 
programs  used  the  five  probe  pressures  and  probe  calibration 
data.  At  the  beginning  of  this  phase  of  the  testing, 
sufficient  probe  calibration  data  were  taken  to  verify  the 
validity  of  the  previous  calibrations.  The  estimated 
accuracy  in  the  calculated  local  angles  of  attack  and 
sideslip,  from  Reference  3,  is  +0.25°. 


Summary  of  Tests 

The  flow-field  survey  tests  are  summarized  in  Table  7. 
Column  1  lists  the  parent  aircraft  configuration  number  from 
Table  2.  The  free-stream  Mach  number,  Mm ,  is  in  column  2 
and  the  probe  angle  of  attack,  ap ,  in  column  3.  The  fourth 
column  indicates  the  Xp  schedule  which  was  used  for  the 
traverse.  These  two  schedules  are  shown  at  the  bottom  of 
each  page  of  the  table.  The  probe  angle  of  attack  listed  in 
the  third  column  is  one  degree  less  than  the  parent  aircraft 
angle  of  attack,  measured  with  respect  to  waterline,  WL, 
zero  shown  in  Figure  1. 

The  remainder  of  Table  7  lists  the  AEDC  run  numbers  for 
the  data  obtained  for  a  particular  parent  aircraft 
configuration,  Mach  number,  probe  angle  of  attack,  and  probe 
position.  The  Xp,  Yp,  Zp  coordinate  system  origin  is  shown 
in  Figure  5.  Blanks  in  the  table  indicate  that  no  data  were 
taken.  The  run  numbers  provide  input  data  for  the  data 
retrieval  computer  program. 
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Use  of  Data  Retrieval  Computer  Program 

The  flow-field  data  are  all  contained  in  one  file  on  a 
magnetic  tape.  A  computer  program  has  been  written  which 
allows  the  flow-field  data  for  run  numbers  listed  in  Table  7 
to  be  retrieved  from  the  tape  and  tabulated.  A  listing  of 
the  computer  program  is  presented  in  Appendix  C.  The  first 
card  of  input  is  in  a  3A4  format  and  contains  the  word 
"FLOW"  beginning  in  column  1.  This  is  used  to  compare  with 
a  header  on  the  tape  to  assure  that  the  correct  tape  is 
mounted.  The  second  input  card  in  an  15  format  specifies 
the  number  of  run  numbers  of  data,  <  100,  to  be  retrieved. 
The  specific  run  numbers  from  Table  7  are  input  on  the  next 
cards,  eight  to  a  card  in  8F10.0  format.  The  run  numbers 
can  be  in  any  order.  The  computer  program  sorts  them  into 
numerical  sequence  before  searching  the  tape. 


Description  of  Tabulated  Data 

There  is  one  page  of  output  for  each  run  number  of  data. 
A  sample  is  shown  in  Figure  12.  The  first  line  of  output 
identifies  the  type  of  data,  in  this  case  flow  field;  the 
AEDC  report  number  documenting  this  test  (Reference  3);  and 
the  ARO  project  number  and  test  number.  The  next  two  lines 
list  the  run  number  and  parent  aircraft  configuration.  The 
following  four  lines  of  output  list  the  free-stream  test 
conditions.  The  next  four  lines  show  the  angles  of  attack 
of  the  parent  aircraft  and  the  probe.  The  remainder  of  the 
page  tabulates  the  probe  position,  in  the  coordinate  system 
of  Figure  5,  and  the  quantities  calculated  from  the  probe 
measurements . 
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FIGURE  12.  Sample  of  Flow-Field  Survey  Data  Tabulation. 
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The  nomenclature  used  in  the  tabulated  output  is 
contained  in  the  following  list. 


ALP  local  angle  of  attack  as  seen  by  the  probe; 

angle  between  the  velocity  component  in  the  Xp, 
Zp  plane  and  that  in  the  Xp  direction,  deg.; 
positive  for  positive  VZ 

ALPHA  PAR  angle  of  attack  of  the  parent  aircraft,  deg. 

ALPHA  PROBE  angle  of  attack  of  the  probe,  deg. 


ALPT  total  angle  of  attack  as  seen  by  the  probe, 

angle  between  probe  longitudinal  axis  and  the 
local  total  velocity  vector,  deg. 


MACH  NO 


fret-stream  Mach  number 


ML 


local  Mach  number  in  the  flow  field 


P 

PHI 


PT 

PTP/PT 


Q 


free-streum  static  pressure,  psfa 

angle  between  the  velocity  component  in  the 
probe  Yp,  Zp  plane  and  the  negative  Zp  axis, 
deg.;  positive  clockwise  viewed  from  the  rear 

free-stream  total  pressure,  psfa 

ratio  of  local  total  pressure  to  free-stream 
total  pressure 

free-stream  dynamic  pressure,  psf 
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OL/Q  ratio  of  local  dynamic  pressure  to  free-stream 

dynamic  pressure 

RE/FT  free-stream  Reynolds  number  per  foot,  ft-'*' 

SIG  local  sidewash  angle  as  seen  by  the  probe; 

angle  between  the  velocity  component  in  the  Xp, 
Yp  plane  and  that  in  the  negative  Xp  direction, 
deg.;  positive  for  positive  VY 

T  free-stream  static  temperature,  °R 

TT  free-stream  total  temperature,  °R 

V  free-stream  velocity,  f t/sec 

VL/V  total  local  velocity  divided  by  the  free-stream 

velocity 

VX/V  local  Xp  velocity  component,  positive  in  the 

negative  Xp  direction,  divided  by  the  freer 
stream  velocity 

VY/V  local  Yp  velocity  component,  positive  in  the 

positive  Yp  direction,  divided  by  the  free- 
stream  velocity 

VZ/V  local  Zp  velocity  component,  positive  in  the 

negative  Zp  direction,  divided  by  the  free- 
stream  velocity 
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XP  position  of  the  probe  static  pressure  orifices 

in  the  Xp  direction  in  the  Xp,  Yp/  Zp  coordi¬ 
nate  system  of  Figure  5,  in. 

YP  position  of  the  probe  longitudinal  axis  in  the 

Yp  direction  in  the  Xp,  Yp,  Zp  coordinate 
system  of  Figure  5,  in. 

ZP  position  of  the  probe  longitudinal  axis  in  the 

Zp  direction  in  the  Xp,  Yp,  Zp  coordinate 
system  of  Figure  5,  in. 


DISCUSSION  OF  ATTACHED  LOADS 


INTRODUCTORY  REMARKS 

A  significant  problem  in  store  loads  is  whether  a  sting- 
supported  store  will  yield  the  attached  store  loads  (forces 
and  moments)  as  it  approaches  the  attached-store  position  as 
close  as  possible  under  practical  testing  conditions  without 
making  contact  with  the  rack.  In  fact,  the  nature  of  the 
interference  forces  on  the  store  for  small  gaps  between  store 
and  rack  is  not  well  understood.  The  purpose  of  this  and  the 
next  section  is  to  elucidate'  the  nature  of  these  loads  as 
revealed  by  the  present  measurements. 

Our  knowledge  of  the  loads  on  a  store  in  close  proximity 
to  the  attached  position  is  not  extensive,  but  light  lias 
been  shed  on  this  subject  by  the  tests  of  Dix  in  Reference 
5.  The  tests  were  of  Mk-Pl  bombs  on  a  r,irt  mounted  on  a 
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model  of  an  F-4C  airplane  and  included  both  attached  loads 
and  Loads  on  a  sting-supported  nonib.  His  results  are  pre¬ 
sented  for  c*s  =  0°  and  do  not  include  parent  aircraft 
component  build-up  results.  For  the  bottom  store,  the 
shoulder  stores  were  present;  for  the  shoulder  store,  the 
opposite  shoulder  store  was  present.  His  data  show  that  the 
captive  loads  cannot  be  obtained  by  extrapolating  CTS  data 
to  zero  Z/D  in  general.  For  some  coefficients  (CLM,  CA, 

CLN ) ,  such  extrapolation  was  better  than  for  others  (CN, 

CY ) .  Large  changes  occurred  for  Z/D  <  0.2  which  could  not 
be  measured  with  the  CTS  system  because  of  limitation  in 
accuracy  of  positioning  the  store. 

In  this  investigation  a  more  extensive  invest igat ion  was 
made  of  the  problem,  and  the  effects  of  angle  of  attack  and 
configuration  build-up  were  measured.  Also  flow-field  data 
in  the  vicinity  of  the  attached  bottom  store  position  were 
taken . 

Before  discussing  the  load  measurements ,  we  will  show  the 
relationship  between  the  attached  and  grid  loads. 

RELATIONSHIP  BETWEEN  ATTACHED  AND  GRID  LOADS 

Stores  S^p  and  SMU  .are  the  ones  tested  both  on 
tii<’  itt  ;ohod- loads  balance  and  the  CTS  sting  supported 
balance.  Data  for  store  SMF  in  combination  with 
conf  inur.it  ion  3  of  Table  2  (F4-C,  inboard  pylon,  TER)  are 

shown  m  Figure  13.  In  Figure  13(a)  the  normal- force 
coefficient  is  shown  versus  7, p/D  for  five  angles  of 
attack  at  Ni  =  0.6.  The  sting  data  extend  down  to 
Zp/D  =  .07.  There  is  a  ] urge  change  in  CN  m  the  range 


59 


NWC  TP  6210 


STORE  VERTICAL  POSITION,  Zp/D 
(a)  Normal-force  coefficient. 


FIGURE  13.  Comparison  of  Attached  and  Grid  Loads  for  Store 
Smf  in  Combination  with  Configuration  3  at  =  0.6. 
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STORE  VERTICAL  POSITION,  Zp/D 
(c)  Pitching-moment  coefficient 


FIGURE  13.  Continued. 


STORE  VERTICAL  POSITION,  Zp/D 
(d)  Yawing-moment  coefficient 
FIGURE  13.  Continued. 
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FIGURE  13.  Concluded 


NWC  TP  6210 


0  <  Zp/D  <  .07  as  indicated  by  the  attached-store  data 
which  cannot  be  reproduced  by  extrapolating  the  sting 
results  to  zero  gap.  This  range  of  Zp/D,  representing 
only  0.05  in.  model  scale,  produces  surprisingly  large 
interference  forces  which  increase  with  angle  of  attack. 

Examination  of  the  other  parts  of  Figure  13  for  the 
other  coefficients  shows  that,  with  the  possible  exception 
of  side-force  coefficient,  extrapolation  of  the  sting  data 
to  Zp  =  0  is  not  accurate.  The  question  naturally 
arises  whether  the  attached  loads  are  accurately  measured. 
Subsequent  data  will  verify  their  repeatability.  It  is  our 
belief  based  on  the  overall  consistency  of  the  data  as  well 
as  the  results  of  Dix  that  the  attached  loads  are  valid. 

The  question  arises  whether  the  above  phenomena  are  also 
valid  at  high  Mach  numbers.  Curves  similar  to  those  of 
Figure  13  are  presented  in  Figure  14  for  store  SMU  in 
combination  with  configuration  3  at  =  0.95.  The  same 
general  behavior  is  seen  as  previously.  In  Figure  15,  the 
normal-force  coefficients  are  shown  for  three  Macn  numoers 
for  store  SMF  in  combination  with  configuration  3. 

These  effects  of  Mach  number  on  the  behavior  for  Zp/D  <  .07 
are  small.  It  appears  that  both  the  sting  and  attached 
load  data  are  repeatable,  and  that  the  strong  interference 
is  not  Mach  number  dependent. 

There  are  adequate  data  to  illustrate  a  number  of 
effects  on  the  way  the  loads  of  the  attached  bottom  store 
vary  with  angle  of  attack.  The  effects  of  the  following 
variables  on  these  attached  store  loads  will  now  be 
illustrated . 


65 


NWC  TP  6210 


0.1  0.2  0.3  0.4 

STORE  VERTICAL  POSITION,  Zp/D 
(c)  Pitching-moment  coefficient 


FIGURE  14.  Continued. 
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(a)  Configuration  (other  than  bottom  store) 

(b)  Mach  number 

(c)  Effect  of  adding  fins  to  shoulder  stores 

(d)  Effect  of  number  of  shoulder  stores 

(e)  Effect  of  bottom  store  configuration 

EFFECT  OF  CONFIGURATION 

Attached  store  loads  (forces  and  moments)  were  measured 
on  four  different  bottom  store  configurations  for 
configurations  3,  4,  5,  6,  and  8.  The  variations  of  loads 
for  store  SAF  at  =  0.6  with  angle  of  attack  are  shown 
in  Figure  16  for  the  above  five  configurations.  Figure 
16(a)  shows  that  at  angles  of  attack  between  -4°  and 
+8°,  configuration  3  yields  the  least  normal  force 
(positive  upward) .  Configurations  4  and  8,  which  only  have 
one  shoulder  store,  yield  more  normal  force;  and 
configurations  5  and  6,  with  two  shoulder  stores,  yield  the 
largest  normal  forces.  What  is  of  particular  interest  is 
that  the  addition  of  fins  to  the  shoulder  stores  has  a  small 
effect  on  normal  force.  Also  at  angles  of  attack  above 
8°,  the  effects  of  configuration  differences  are  much 
less  than  at  lower  angles  of  attack. 

The  side- force  variation  with  ag  shown  in  Figure  16(b) 
shows  that  the  side  force  is  relatively  insensitive  to 
configuration  especially  for  angles  of  attack  less  than  8°. 
Positive  side  force  is  directed  towards  the  fuselage  from 
the  left-wing  inboard  pylon. 

Pitching-moment  coefficients  shown  in  Figure  16(c)  are 
generally  nose  down.  For  angles  of  attack  less  than  about 
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(c)  Pitching-moment  coefficient. 
FIGURE  16.  Continued. 
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(e)  Rolling-moment  coefficient 
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6°,  the  effect  of  configuration  on  the  magnitude  of  the 
pitching-moment  coefficient  is  similar  to  its  effect  on 
normal  force.  However  in  the  range  above  6°  to  8°, 
configurations  5  and  6  have  smaller  pitching  moments  than 
configurations  4  and  8.  Again  adding  fins  to  the  shoulder 
stores  causes  negligible  effect. 

The  yawing-moment  coefficient  is  slightly  nonlinear  with 
from  -4°  to  16°  and  is  positive.  There  are  significant 
differences  due  to  configuration  changes,  but  the  addition 
of  fins  to  the  shoulder  store  has  a  negligible  effect. 

The  rolling-moment  coefficients  in  Figure  16(e)  are 
positive,  nonlinear,  and  have  maxima  around  6°  to  8°.  The 
configuration  effects  are  significant.  The  addition  of  the 
fins  to  the  shoulder  stores  causes  a  significant  effect  in 
the  higher  angle  range. 

The  same  kind  of  data  as  Figure  16  for  M^  =  0.95  show 
the  same  general  qualitative  effects. 

EFFECT  OF  MACH  NUMBER 

A  set  of  curves  for  the  loads  on  SAF  for 
configuration  6  at  M  =  0.6,  0.7,  0.8,  0.9  and  0.95  in 
Figure  17  shows  that  Mach  number  has  its  largest  effect  on 
normal  force.  Mach  number  effects  on  the  rest  of  the 
coefficients  are  not  large. 
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FIGURE  17.  Effect  of  Mach  Number  on  the  Attached 
Loads  of  Bottom  Store  Saf  for  Configuration  6. 
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FIGURE  17 
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(c)  Pitching-moment  coefficient. 


FIGURE  17.  Continued 


NWC  TP  6210 


1  f — V;  \ 


nD  ‘jiNaiDi Jiiaoo  iNawow-onmoa 


83 


STORE  ANGLE-OF- ATTACK,  ag,  DEG 
(e)  Rolling-moment  coefficient. 


NWC  TP  6210 


EFFECT  OF  ADDING  FINS  TO  SHOULDER  STORES 

The  effect  of  adding  fins  to  the  shoulder  stores  on  the 
bottom  store  loads  is  seen  by  comparing  results  for 
configurations  4  and  8  or  configurations  5  and  6  (latter 
case  corresponding  to  flight).  The  loads  on  store  SAF 
are  shown  in  Figure  18  for  Mm  =  0.6  for  configurations  5  and 
6.  The  data  show  negligible  effect  of  adding  the  fins  to 
the  stores  except  for  rolling-moment  coefficient  and  for 
pitching-moment  coefficient  at  high  angle  of  attack.  At 
=  0.95  the  same  kind  of  results  show  negligible  effect  of 
adding  fins  to  the  shoulder  stores  except  on  CN  and  CLN  at 
negative  angles  of  attack,  and  on  rolling  moment. 


EFFECT  OF  NUMBER  OF  SHOULDER  STORES 

By  comparing  store  loads  for  configurations  4  and  5  or  6 
and  8,  the  effect  of  the  number  of  shoulder  stores  on  the 
bottom  store  loads  can  be  determined.  Figure  19  presents 
the  loads  on  SAF  at  =  0.6  for  configurations  4  and  5. 

The  addition  of  a  second  shoulder  store  has  the  effect  of 
increasing  normal-force  coefficient  at  low  angles  of  attack 
and  decreasing  it  at  high  angles  of  attack  with 
corresponding  changes  in  the  pi tch ing-mornont  coefficient. 
While  the  addition  of  the  second  shoulder  store  has  a  small 
effect  on  side-force  coefficient  except  at  negative  angle  of 
attack,  it  does  have  a  significant  effect  on  yawing-moment 
coefficient.  Rolling-moment  coefficient  is  also  changed. 

At  =  0.95  the  data  exhibit  the  same  qualitative 
effects  as  at  fl  =  0.6. 
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FIGURE  18.  Effect  of  Adding  Fins  to  Shoulder  Stores  on  Attached 
Loads  of  Bottom  Store  Sap  at  =  0.6. 
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STORE  ANGLE-OF- ATTACK,  «s,  DEG 
(b)  Side-force  coefficient. 
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FIGURE  19.  Effect  of  Number  of  Shoulder  Stores  on  the  Attached 
Loads  of  the  Bottom  Store  Sid  at  M  =  0.6. 
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STORE  ANGLE-OF- ATTACK,  ag,  DEG 
(c)  Pitching-moment  coefficient. 
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EFFECT  OF  BOTTOM  STORE  CONFIGURATION 

It  is  of  interest  to  know  how  the  four  stores  tested, 
SAF'  SAU'  SMF'  and  SMU'  varY  i-n  their  attached  loads.  Such 
results  for  =0.6  are  shown  in  Figure  20  for  the  stores 
attached  to  configuration  6.  The  stores  are  all  different 
with  respect  to  normal-force  and  pitching-moment 
characteristics.  With  regard  to  side-force  and 
yawing-moment  coef f icients ,  the  finned  stores  are  nearly  the 
same  and  the  unfinned  stores  are  also  nearly  the  same. 
Significant  difference  exist  between  the  finned  stores  with 
regard  to  rolling-moment  coefficient.  The  measured  rolling 
moments  for  the  unfinned  stores  are  an  indication  of  the 
accuracy  of  the  measurements. 


DISCUSSION  OF  GRID  LOADS 

CONFIGURATION  EFFECTS  AT  ug  =  0° 

The  grid  loads  on  the  bottom  store  as  it  moves  downward 
from  close  proximity  to  the  rack  are  influenced  by  angle  of 
attack,  configuration,  M.icli  number,  addition  of  fins  to  the 
shoulder  stores,  and  the  number  of  shoulder  stores.  Figure 
21  shows  the  loads  on  store  at  M^  =  0.6  as 

influenced  by  configuration  and  store  vertical  position. 
Only  the  first  1.25  diameter  of  store  vertical  position  are 
shown . 

Examination  of  the  normal - force  results  shows  the  least 
normal  force  for  configuration  2  with  a  small  positive 
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FIGURE  20.  Dependence  of  Attached  Loads  of  Bottom  Store  on 
Its  Configuration  for  Airplane  Configurations  at  =  0.6. 
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(c)  Pitching-moment  coefficient. 
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STORE  VERTICAL  POSITION,  Zp/D 
a)  Normal-force  coefficient. 

FIGURE  21.  Grid  Loads  on  Store  Sj^p  as  Influenced  by 
Airplane  Configuration;  MTO  -  0.6,  as  =  0°. 
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(b)  Side-force 
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STORE  VERTICAL  POSITION,  Zp/D 
(c)  Pitching-raoment  coefficient. 
FIGURE  21.  Continued. 
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increment  due  to  addition  of  the  TER.  Addition  of  one 
shoulder  store  contributes  positive  increments  about  twice 
those  due  to  the  TER.  Addition  of  the  other  shoulder  store 
adds  an  even  larger  positive  increment  to  the  normal  force. 
Addition  of  the  fins  to  the  shoulder  stores  adds  further 
small  positive  increments  in  normal  force.  We  conclude  that 
the  addition  of  the  shoulder  stores  are  most  important, 
followed  by  the  addition  of  the  TER,  and  finally  the 
addition  of  the  fins  to  the  shoulder  stores. 

The  side-force  results  in  Figure  21(b)  show  that  the 
largest  effects  of  configuration  occur  at  small  Zp/D. 

The  effect  of  adding  the  second  shoulder  store  is  to 
almost  cancel  the  increment  due  to  adding  the  first. 

The  effect  of  adding  fins  to  the  shoulder  stores  is  again 
fairly  smal 1 . 

The  pitching-moment  results  show  the  same  qualitative 

effects  as  the  normal-force  results.  The  yawing-moment 

results  show  large  configuration  effects  at  small  Zp/D 

(on  an  expanded  scale)  and  small  effects  at  large  Zp/D. 

The  store  alone  has  a  rolling-moment  coefficient  of  .037  at 

«  =  n°  because  of  fin  cant.  Additional  rolling  moments  of 

about  the  same  magnitude  can  be  developed  at  small  Zp/D 

because  of  configurational  differences.  At  these  small 

values  of  Zp/D  increments  in  CLL  as  large  as  0.01  can  be 

developed  by  the  addition  of  fins  to  the  shoulder  stores. 

The  store  alone  has  rolling-moment  coefficients  in  the  range 

of  .026-. 037  at  <f>  =0°  and  .026-. 040  at*  =  22.5° 

s  s 

for  0°  <  ag  <  8°  principally  by  virtue  of  the  cant 
of  the  tail  fins.  (Note  that  $  =  0°  corresponds  to  the 

X  tail  roll  orientation.) 
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CONFIGURATION  EFFECTS  AT  a  =8° 

s 

Figure  22  shows  the  same  kind  of  data  as  in  Figure  21 

except  for  a  =  8°  instead  of  a  =  0°.  The  results 
s  s 

of  Figure  22(a)  show  that  the  normal-force  coefficient 

increases  as  the  bottom  store  moves  away  from  the  rack 

rather  than  remains  constant  or  decreases  as  at  a  =  0°. 

s 

Also  the  addition  of  the  outboard  shoulder  store  to 
configuration  3  now  reduces  the  normal  force  rather  than 
increases  it.  The  effect  of  adding  the  fins  to  the  shoulder 
store  affects  the  bottom  store  normal  force  only  for  small 
values  of  Zp/D. 

The  side-force  coefficients  of  Figure  22(b)  are  much 

higher  than  those  of  Figure  22(a),  and  the  effects  of 

configuration  changes  are  generally  much  less.  There  is  not 

much  qualitative  difference  in  the  pitching-moment 

coefficient  behavior  at  =  8°,  Figure  22(c),  from  that 

at  =  0°.  For  the  yawing-moment  coefficient,  lower 

values  prevail  at  a  =  8°,  Figure  22(d),  than  at  a  = 

s  s 

0°,  with  some  qualitative  differences  in  behavior.  The 

peak  rolling-moment  coefficients  at  =  8°  do  not 

exceed  those  at  u  =  0°  by  much. 

s 

Vlith  regard  to  the  effect  of  adding  fins  to  the  shoulder 
stores  on  the  results,  compare  configurations  5  and  0  in 
Figure  22.  Adding  the  fins  generally  causes  a  sin  a  1  1  but 
significant  effect  at  small  values  of  Zp/'D.  The  effect 
is  usually  much  smaller  than  that  due  to  adding  t ue  :dw  wider 
stores  to  the  rack. 
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EFFECT  OF  MACH  NUMBER 

The  effects  of  Mach  number  on  the  grid  loads  of  bottom 

store  SMF  are  shown  in  Figure  23  for  configuration  6  at 

a  =  0°.  The  effects  of  Mach  number  on  normal-force 
s 

coefficient  are  not  large  for  M^  =  0.6-0. 8,  for  the  entire 
range  of  Zp/D.  At  low  values  of  this  parameter,  the  spread 
between  the  M  =0.6  and  0.95  data  becomes  fairly  significant 

oo 

This  result  also  is  true  of  the  attached  load  data  seen  in 
Figure  17(a)  for  the  same  test  conditions  but  for  store  SAF. 

The  side-force  results  of  Figure  23(b)  show  small 
effects  of  Mach  number  for  all  values  of  Zp/D,  when  the 
expanded  scale  on  the  figure  is  taken  into  consideration. 

The  pitching-moment  results  of  Figure  23(c)  show  a 
systematic  effect  of  Mach  number.  The  data  on  yawing  moment 
show  little  effect  of  Mach  number  at  high  values  of 
Zp/D,  but  at  low  values  the  data  for  M^  =  0.6,  0.7,  and 
0.8  form  one  trend  and  those  for  0.9  and  0.95  another  trend. 
The  rolling  moment  data  all  show  the  same  trend  with 
Zp/D  with  secondary  variations  due  to  M^ . 

It  is  of  interest  to  contrast  the  effect  of  Mach  number 

on  the  loads  at  a  =  8°  shown  in  Figure  24  with  those  at 

s 

a  =  0°  shown  in  Figure  23.  The  effects  of  M  at  a  = 

S  °o  S 

8°  on  normal-force  coefficient  are  about  the  same 

magnitude  as  they  are  at  ag  =  0°.  The  effects  of  Mach 

number  on  side-force  coefficient  at  a  =8°  are 

s 

proportionally  about  the  same  as  at  ocs  =  0°,  but  the 
increments  are  much  larger  since  the  side-force  coefficients 
are  much  larger. 
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STORE  VERTICAL  POSITION,  Zp/D 

(a)  Normal-force  coefficient. 

FIGURE  23.  Grid  Loads  on  Store  as  Influenced 

by  Mach  Number;  Configuration  6,  a  =  0°. 
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STORE  VERTICAL  POSITION,  Zp/D 
(c)  Pitching-moment  coefficient. 
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SYM  MACH  RUN 
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STORE  VERTICAL  POSITION,  Zp/D 
(c)  Pitching-moment  coefficien 
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A  systematic  effect  of  on  the  pitching-moment 

coefficient  is  manifest  at  a  =  8°  just  as  at  a  = 

0°  but  the  qualitative  effect  of  Zp/D  is  different 

in  the  two  cases.  The  yawing-moment  qualitative  variation 

with  Zp/D  is  quite  different  at  ag  =  8°  than  at  = 

0°.  The  effects  of  M  on  the  data  are  least  for  Mm  = 

0.6-0. 8  and  larger  for  M  =  0.8-0.95.  As  at  a  =  0°, 

the  trend  of  rolling-moment  coefficient  with  Zp/D  at 

a  =8°  dominates  the  behavior  with  a  secondary  effect  of 
s 

superimposed  on  the  Zp/D  trend. 

EFFECT  OF  ADDING  FINS  TO  SHOULDER  STORES 

It  is  a  complication  to  model  the  shoulder  store  fins  in 
the  computer  program.  It  is  therefore  of  interest  to  see 
what  is  lost  in  accuracy  by  not  doing  so.  A  measure  of  the 
inaccuracies  incurred  can  be  obtained  by  comparing  the  loads 
on  the  bottom  store  at  identical  positions  in  configurations 
5  and  6.  These  data  for  SMp  at  =  0°  are  shown  in 
Fiqure  25  for  M  =0.6.  The  fin  effects  on  normal-force 
coefficient  are  a  maximum  of  about  10  percent  at  small 
Zp/D.  The  maximum  fin  effect  on  side  force  coefficient 
is  about  .02,  also  at  very  low  Zp/D.  There  is  about  a 
maximum  10  percent  effect  on  pitching-moment  coefficient. 

The  yawing-moment  coefficients  are  quite  small  being  both 
negative  and  positive.  An  error  of  about  0.1  in  the 
coefficient  occurs  at  very  low  values  of  Zp/D.  At  low 
Zp/D  the  rolling-moment  coefficient  difference  is  about 
0.01. 


It  is  of  interest  to  look  at  other  cases  to  see  the 
effect  of  adding  the  fins  to  the  shoulder  stores  as 
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FIGURE  25.  Grid  Loads  on  Store  smf  as  Influenced  by 
Addition  of  Fins  to  Shoulder  Stores;  14  =0.6,  a  =  0°. 
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STORE  VERTICAL  POSITION 
(e)  Rolling-moment  coeff 
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influenced  by  and  Qig.  The  following  table  gives  some 
results  from  figures  not  shown. 


M  =0.6 

00 

a  =  8° 
s 

CN: 

.05  in  .20 

CY: 

negligible 

CLM : 

0.15  in  0.4 

CLN: 

10% 

CLL : 

negligible 

M  =0.95 

p 

ii 

o 

o 

00 

s 

CN: 

10% 

CY: 

.015  in  0.11 

CLM: 

0.25  in  1.75 

CLN: 

negligible 

CLL: 

.008  in  .06 

=  0.95 

> 

0 

CO 

II 

Ul 

3 

CN: 

negligible 

CY: 

negligible 

CLM: 

0 . 2  out  of  0.8 

CLN: 

.08  out  of  0.7 

CLL: 

negl igible 

Since  the  effect  of  the  shoulder  store  fins  is  usually 
greatest  for  small  Zp/D,  it  is  probable  that  in  a  store 
separation  trajectory  the  effect  of  those  fins  can  be 
neglected  for  many  cases.  However  for  accurate  attached 
load  calculations  they  should  probably  be  included. 
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EFFECT  OF  SHOULDER  STORES 

Figure  26  compares  the  loads  on  store  SMF  for 
configurations  3,  5  and  6.  These  results  show  the 

importance  of  the  shoulder  stores  on  the  loads  of  the  bottom 
store  and  how  far  out  the  shoulder  stores  affect  the  loads. 
For  normal-force  and  pitching-moment  coefficients  the  second 
store  adds  an  increment  about  as  large  as  that  due  to 
addition  of  one  store,  and  the  interference  effects  are 
significant  out  to  4-6  store  diameters.  Comparable  results 
are  shown  for  side  force  and  yawing  moment  except  that  the 
interference  effects  are  only  significant  to  1-2  store 
diameters.  The  rolling  moment  results  indicate  large 
effects  only  for  small  Zp/D  values,  and  show  a  difference  of 
magnitude  of  .005  for  Zp/D  >  1.0. 

The  range  of  the  interference  effects  of  the  shoulder 
stores  on  the  bottom  store  is  judged  to  be  large  enough  that 
it  can  have  significant  effects  on  trajectory  calculations. 
Such  stores  also  must  be  accounted  for  in  determining 
attached  loads. 


EFFECT  OF  SWAY  BRACES 

Measurements  have  been  made  to  determine  the  importance 

on  store  grid  loads  of  sway  braces  such  as  shown  in  Figure 

5.  Measurements,  with  and  without  sway  braces,  are 

indicated  in  Table  6.  Comparing  the  grid  data  for  the 

bottom  store  of  configuration  6  at  M  =  0.6,  0.8,  0.9,  and 

0.95  at  a  =  0°  with  and  without  sway  braces,  we  find 
s 

their  effects  are  generally  negligible  except  in  a  few  cases 
for  close  proximity  to  the  attached  position  for  CN  and  CLM . 
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FIGURE  26.  Grid  Loads  on  Store  smf  as  Influenced 
by  Number  of  Shoulder  Stores;  M.^  =  0.6,  as  =  0  . 
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(e)  Rolling-moment  coefficient. 
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For  CN  of  0.6,  an  effect  of  0.1  for  Zp/D  <  0.1  is  about 
the  maximum.  For  a  CM  of  1.2  a  change  of  0.1  to  0 . 2  is  a 
maximum.  These  changes  attenuate  to  very  small  values  by 
Zp/D  =  0.5.  These  results  are  representative  of  all 
Mach  numbers. 

Dix  in  Reference  6  examined  the  effects  of  sway 

braces  on  the  attached  loads  for  the  same  configuration. 

At  a  =  0°  his  results  were  the  same  as  described  above 
s 

except  for  CM.  He  found  increments  in  CM  as  high  as 

0.6. 


DISCUSSION  OF  FLOW  FIELD  MEASUREMENTS 

The  nature  of  the  flow  field  is  influenced  by  the  same 
parameters  which  affect  the  store  loads.  They  include 
vertical  position,  angle  of  attack,  configuration,  Mach 
number,  number  of  shoulder  stores,  and  addition  of  fins  to 
the  shoulder  stores.  The  quantities  shown  in  the  ensuing 
figures  are  VZ/V  and  VY/V,  which  are,  respectively,  the 
upwash  and  sidewash  velocities  divided  by  the  free-stream 
velocity.  For  the  purposes  of  this  report,  these  quantities 
can  be  equated  to  the  upwash  and  sidewash  angles  in  radians 
within  a  few  percent  in  all  cases.  The  angles  are  measured 
in  the  attached  store  or  probe  axis  system  shown  in 
Figure  5. 
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EFFECT  OF  VERTICAL  POSITION  AND  ANGLE  OF  ATTACK 


Figure  27  shows  measured  upwash  and  sidewash  under  the 
rack  of  configuration  3  at  =  0.6  at  Zp/D  of  0  and  1.0 
for  both  0°  and  16°  angles  of  attack.  It  is  noted  that 
the  change  in  upwash  between  these  vertical  positions  for 
both  angles  of  attack  is  not  large.  The  sidewash  angle  in 
Figure  27(b)  tends  to  show  more  variation  between  positions. 

The  results  which  follow  are  for  Zp/D  =  0. 


EFFECT  OF  CONFIGURATION 

The  effects  of  configuration  will  first  be  examined  with 
regards  to  adding  the  pylon  and  then  the  rack  to  the  clean 
airplane  by  comparing  data  for  conf igurations  1,  2,  and  3. 
The  upwash  and  sidewash  for  =  0.6  and  ag  =  0°  are  shown 
in  Figure  28  for  these  configurations  as  well  as  others. 
Adding  the  pylon  to  the  clean  airplane  hardly  changes  the 
upwash,  but  adding  the  rack  causes  a  significant  effect; 
as  much  as  2°.  The  effect  of  adding  the  pylon  on  sidewash 
has  a  maximum  effect  of  about  0.5°  while  adding  the  rack 
has  about  three  times  this  effect. 

Examining  similar  results  in  Figure  29  for  ag  =  4°,  we 
note  results  of  about  the  same  magnitude  for  upwash  angle. 
The  results  for  sidewash  in  Figure  29(b)  are  similar  to 
those  for  =  0°  but  are  of  larger  magnitude. 

It  seems  desirable  to  model  the  pylon  and  necessary  to 
model  the  rack. 
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The  effects  of  adding  shoulder  stores  to  the  rack  and 
pylon  and  fins  to  the  shoulder  stores  can  be  seen  by 
comparing  configurations  3,  5,  and  6  in  Figures  28  and  29. 

It  is  noted  that  configurations  5  and  6  cause  large  changes 
in  upwash  from  configuration  3  in  Figure  28(a)  but  they  have 
virtually  identical  upwash  for  0  <  Xp  <  -5.  The  Xp  =  -5 
position  is  near  the  fin  leading  edge,  and  only  behind 
this  point  is  the  upwash  field  due  to  the  addition  of  the 
fins  of  significance. 

Adding  the  shoulder  stores  to  the  rack  causes 
significant  increments  in  sidewash  but  these  are  usually 
less  than  those  caused  by  adding  the  rack  to  the  pylon.  The 
further  addition  of  the  fins  causes  small  sidewash 
increments  which  are  felt  in  front  of  and  behind  the  fins. 

At  c*s  =  4°,  the  same  general  conclusion  holds  for 
upwash,  while  the  maximum  increment  in  sidewash  is  about 
1.0°  due  to  adding  shoulder  stores  and  another  0.5° 
due  to  addition  of  fins. 

It  thus  appears  that  addition  of  the  rack  has 
significant  effects  on  upwash  and  sidewash.  Addition  of  a 
pair  of  shoulder  stores  also  has  significant  effects  on  the 
upwash  and  sidewash.  Addition  of  fins  to  the  shoulder 
stores  has  little  effect  on  upwash,  and  then  only  behind  the 
fins.  It  has  a  small  effect  on  sidewash  in  front  of  and 
behind  the  fins. 

EFFECT  OF  MACH  NUMBER 

The  effect  of  Mach  number  on  the  upwash  and  sidewash  at 
=  4°  is  shown  for  configuration  3  in  Figure  30.  There 
is  an  effect  of  Mro  on  upwash  in  the  region  of  large  axial 
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FIGURE  30.  Effect  of  Mach  Number  on  Flow  Angles  at 
Zn/D  =  0  for  aa  =  4°;  Configuration  3. 
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gradients,  but  there  is  a  small  effect  on  sidewash  for  the 
range  0-6  <  Me  _<  0.95.  Similar  data  are  presented  for 
configuration  6  in  Figure  31.  There  are  small  effects  of 
on  upwash  and  sidewash  up  to  =  O.S.  There  is  a 
maximum  change  in  sidewash  angle  of  about  1°  over  a  small 
Xp  range  due  to  change  in  from  0.8  to  0.95. 

EFFECT  OF  NUMBER  OF  SHOULDER  STORES 


The  effect  of  adding  the  outboard  shoulder  store 
(including  fins)  to  the  rack  is  found  by  comparing  config¬ 
urations  3  and  8,  and  the  effect  of  adding  both  shoulder 
stores  (including  fins)  to  the  rack  is  found  by  comparing 
configurations  3  and  6.  These  effects  are  shown  on  upwash 
and  sidewash  for  Zp/D  =  0  at  =  0.6  in  Figure  32  for  = 

0°  and  in  Figure  33  for  a  =  4°.  Examining  Figure  32(a), 

s 

we  can  see  that  adding  both  stores  causes  generally  more 
than  twice  the  effect  on  upwash  of  adding  one  store.  The 
sidewash  angles,  Figure  32(b),  are  negative,  and  for  the 
present  measurements  under  the  left  wing,  the  flow  is 
directed  outboard.  The  effect  of  adding  two  shoulder  stores 
is  not  to  change  the  sidewash  angle  for  the  first  two  inches 
(one-third  of  the  store  length),  while  one  shoulder  store 
decreases  the  sidewash  angle  substantially.  There  is  a 
clearcut  difference  between  the  effects  of  adding  one  or  two 
shoulder  stores  on  the  sidewash. 

The  upwash  results  for  ag  =  4°  in  Figure  33(a)  show  that 
adding  the  outboard  shoulder  store  has  a  greater  effect  over 
the  forward  half  of  the  survey,  and  that  adding  the  inboard 
store  has  the  greater  effect  over  the  rear  half  of  the  survey. 
With  regard  to  sidewash,  the  results  noted  in  Figure  32(a) 
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FIGURE  31.  Concluded. 
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PROBE  AXIAL  POSITION,  Xp,  IN. 

(a)  Upwash  angle. 

FIGURE  33.  Effects  of  Adding  Shoulder  Stores  to  Rack  and  of  Adding  Fins 
to  Shoulder  Stores  on  Flow  Angles  at  Zp/D  =  0  and  ^  =  0.6;  as  =  4°. 
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PROBE  AXIAL  POSITION,  Xp,  IN. 
(b)  Sidewash  angle. 
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for  0  >  Xp  >  -2  are  repeated,  but  the  effect  of  adding  either 
or  both  stores  on  sidewash  is  small  for  Xp  <  -3. 


EFFECT  OF  ADDING  FINS  TO  SHOULDER  STORES 


Comparison  of  configurations  5  and  6  shows  the  effect  of 
adding  fins  to  the  shoulder  stores.  These  results  are  given 
for  Mm  =  0.6  in  Figure  34  for  a  =  0°  and  in  Figure  35 
for  as  =  4°.  With  respect  to  upwash  there  are  no  large 
effects  of  adding  the  fins.  There  is  about  a  1°  change 
in  upwash  at  ag  =  0°  behind  the  fins  of  the  store  due  to 
the  addition  of  the  fins.  The  changes  in  sidewash  are 
small,  not  exceeding  1/4  to  1/2  degree.  At  ag  =  4°, 
there  is  little  change  in  upwash  and  about  a  maximum  of 
0.4°  in  sidewash. 


COMPARISON  BETWEEN  EXPERIMENT  AND  THEORY 


GENERAL  APPROACH 

The  basic  objective  of  the  experimental  investigation  is 
to  provide  data  for  validating  the  computer  program  of 
Reference  1  and  for  providing  insight  into  methods  for 
upgrading  the  computer  program.  Accordingly  the  comparisons 
between  theory  and  experiment  are  directed  toward  these 
ends.  In  order  to  make  such  comparisons  meaningful,  it  is 
of  interest  to  describe  the  main  assumptions  that  are  made 
in  the  computer  program  especially  with  respect  to  pylon, 
rack,  and  stores  and  to  suggest  areas  where  refinements  in 
the  computer  program  may  appear  necessary. 
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FIGURE  34.  Effect  of  Adding  Fins  to  Shoulder  Stores  on 
Flow  Angles  at  Zp/D  =  0  for  =  0.6;  a  =  0°. 
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ASSUMPTIONS  OF  THE  COMPUTER  PROGRAM 

The  computer  program  is  based  on  linear  theory, 

Laplace's  equation  for  the  flow  field  with  the 
Glauert-Prandtl  theory  to  account  for  compressibility.  The 
fuselage,  which  can  be  noncircular,  is  modeled  by  using 
sources,  doublets,  and  high-order  solutions  along  the  axis 
of  the  fuselage  to  satisfy  its  boundary  conditions,  with  a 
vortex- lattice  layout  on  the  wing  to  satisfy  the  wing 
boundary  conditions.  The  vortex-lattice  system  is  imaged  in 
the  fuselage  (assumed  circular).  Wing  thickness  is  modeled 
by  source  panels.  Account  is  taken  of  airflow  through 
inlets  and  ducts  by  changing  the  effective  cross-sectional 
area  distribution  of  the  body. 

The  pylon  is  modeled  with  regard  to  thickness  by  source 
panels  and  with  regard  to  the  normal  velocity  boundary 
condition  by  a  vortex-lattice  layout. 

The  rack  is  modeled  for  volume  by  a  body  of  revolution. 

The  stores  are  modeled  for  volume  by  three-dimensional 
source  distributions  along  their  axes.  No  mutual  inter¬ 
ference  between  stores  or  between  stores  and  rack  is 
accounted  for.  No  doublets  to  model  store  angle  of  attack 
distributions  are  included.  The  tail  fins  of  the  shoulder 
stores  are  not  modeled. 

The  forces  on  the  ejected  store  are  calculated  by 
slender-body  theory.  Both  upwash  and  sidewash  distribu¬ 
tions  are  taken  into  account.  Also  a  loading  due  to 
buoyancy  is  included  in  the  calculation.  If  the  flow 
separates  at  some  axial  station,  crossflow  drag  theory  is 
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used  to  calculate  the  loading  downstream  of  separation. 

The  tail  fin  contributions  are  determined  using  the  spanwise 
variation  of  induced  downwash  across  the  tail  span  together 
with  reverse-flow  theorems  in  a  method  which  has  accuracy 
nearly  equivalent  to  full  linear  theory. 

POSSIBLE  SHORTCOMINGS  OF  COMPUTER  PROGRAM 

In  the  ensuing  remarks  a  distinction  will  be  made 
between  the  problem  of  determining  the  loads  on  a  separating 
store  sufficiently  accurately  to  calculate  its  separation 
trajectory  and  the  problem  of  accurately  predicting  attached 
loads.  The  large  changes  in  loads  in  the  first  few  tenths 
of  diameter  of  store  downward  travel  do  not  need  to  be 
accurately  predicted  for  the  first  problem  because  the  store 
spends  very  little  of  its  total  time  under  the  influence  of 
the  aircraft  in  this  region.  With  regard  to  the  second 
problem  it  should  also  be  borne  in  mind  that  the  attached 
aerodynamic  loads  are  also  augmented  by  inertial  loads,  and 
therefore  represent  only  part  of  the  maximum  loads. 

The  questions  that  need  examination  with  respect  to  the 
computer  program  are  largely  those  of  whether  the  airplane, 
pylon,  rack,  and  stores  are  individually  adequately  modeled, 
and  whether  their  mutual  interferences  are  properly  accounted 
for.  Additional  questions  arise  concerning  the  limits  of 
the  computer  program  because  of  nonlinear  effects  of  angle 
of  attack  and  Mach  number.  Some  remarks  concerning  the 
pylon- rack-store  interference  problem  should  be  helpful  in 
the  subsequent  comparisons  between  experiment  and  theory. 
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With  regards  to  Figure  5,  the  pylon  is  a  lifting  surface 
which  should  be  modeled  fairly  accurately  by  vortex-lattice 
methods.  Its  effect  on  the  store  is  minimized  by  the  fact 
that  the  store  is  fairly  far  away  because  of  the  intervening 
rack.  At  angle  of  attack  =  4°,  Figure  29(b)  shows  about 
3°  of  sidewash  at  the  store  position.  At  the  rack,  the 
sidewash  would  be  still  higher.  The  rack  is  modeled  as  a 
body  of  revolution,  and  no  doublets  to  account  for  downwash 
or  sidewash  at  the  rack  are  included  in  the  model. 

Figure  3  shows  a  detailed  model  of  the  TER.  The  first 
0.8  in.  of  the  rack  has  cross  sections  similar  to  the  upper 
half  of  the  cross-sectional  shape  shown.  This  part  of  the 
rack  is  probably  modeled  sufficiently  accurately  by  a  body 
of  revolution.  From  0.8  in.  to  about  2.2  in.,  the  rack  has 
the  A-A  cross  section  shown.  In  effect  a  secondary  pylon 
exists  beneath  the  rack  which  is  in  close  proximity  to  the 
store.  This  secondary  pylon  is  not  properly  accounted  for 
by  modeling  the  rack  as  a  body  of  revolution.  The  rear  part 
of  the  rack  is  hexagonal,  and  is  well  modeled  by  a  body  of 
revolution.  The  rack  as  a  body  of  revolution  is  over  a  rack 
diameter  away  from  the  lower  store.  Since  the  flow  field  of 
three-dimensional  sources  drops  off  with  radius  as  r 
and  that  due  to  doublets  as  r  ,  it  is  to  be  expected 
that  the  source  terms  due  to  the  rack  will  dominate  the 
lower  store.  For  the  shoulder  stores  it  may  not  be  true 
that  this  is  the  case  because  of  the  close  proximity  of  the 
rack  to  the  shoulder  stores. 

Consider  now  the  shoulder  store  modeling.  In  the 
computer  codes  they  are  modeled  only  by  source  distributions 
with  no  mutual  interference  between  each  other  or  between 
rack  and  lower  store.  There  are  upwash  and  sidewash  angles 
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at  the  shoulder  store  locations  which  should  be  modeled  for 
attached  loads.  Also  the  mutual  interference  among  all 
stores  and  the  rack  should  be  accounted  for. 

In  determining  the  loads  on  the  lower  store,  the  flow 
field  induced  directly  by  the  rack  and  shoulder  stores  is 
taken  into  account  with  no  mutual  interference  between  these 
components.  Such  mutual  interference  should  probably  be 
accounted  for  to  obtain  better  estimates  of  the  induced  flow 
field  at  the  bottom  store  position.  When  the  bottom  store 
is  introduced,  it  produces  additional  mutual  inte-  ference 
between  itself  and  the  rack  and  shoulder  stores.  This 
mutual  interference  arises  from  bottom  store  source  and 
doublet  distributions.  It  seems,  that  as  a  minimum,  mutual 
interference  among  the  sources  should  be  accounted  for,  and 
possibly  also  doublets  near  the  attached  position. 

Finally,  tail  fins  on  the  shoulder  stores  are  not 
accounted  for  in  the  computer  program.  Modeling  them  is 
probably  of  secondary  importance  compared  to  the  above 
shortcomings  based  on  the  small  differences  seen  between 
configurations  5  and  6. 

COMPARISON  BETWEEN  DATA  AND  THEORY 
FOR  CLEAN  AIRPLANE 

In  the  ensuing  comparisons,  the  effects  of  adding  the 
pylon  to  the  clean  airplane,  the  effect  of  adding  the  rack 
to  the  pylon,  and  the  effect  of  adding  the  shoulder  stores 
to  the  rack  will  be  isolated  and  compared  with  theory. 

These  increments  are  all  to  be  added  to  the  clean  airplane 
characteristics  as  a  base  configuration.  It  is  therefore 
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of  interest  to  examine  the  comparison  between  experiment  and 
theory  for  the  clean  airplane. 

In  Figure  36  the  downwash  and  sidewash  angles  are  shown 
along  the  Zp/D  =  0  location  for  =  0.6  for  various 
angles  of  attack.  It  is  noted  that  the  experimental  upwash 
angle  is  predicted  fairly  well  up  to  about  8°  angle  of 
attack  although  it  is  consistently  less  than  theory. 
Significantly  large  deviations  between  experiment  and  theory 
occur  by  the  time  a  =  16°  is  reached.  The  sidewash  is 
predicted  well  up  to  a =  8°.  The  deviations  that  do  exist 
at  the  higher  angles  could  be  due  to  a  number  of 
considerations  such  as  departure  of  the  aircraft  lift  from 
the  range  of  linearity,  imprecise  model  dimensions,  tunnel 
flow  angularity,  and  mismatch  of  the  engine  mass  flow  ratio 
between  the  theoretical  model  and  the  wind-tunnel  model . 

Before  comparing  the  theory  and  experiment  for  the 
forces  and  moments  measured  under  the  conditions  of  Figure 
36,  it  is  of  interest  to  describe  the  characteristics  of 
store  SMF  used  in  the  comparison.  Test  data  for  the 
store  with  and  without  tail  fins  in  the  X  configuration  were 
used  to  determine  the  increment  in  normal  force  due  to 
adding  the  tail  fins  to  the  body.  The  increment  was  divided 
between  tail  fins  and  the  body  in  accordance  with  the  Kw 
and  Kg  interference  factors  described  in  Reference  7. 

The  body  alone  has  a  boattail  with  possible  separation.  A 
crossflow  drag  coefficient  was  selected  for  the  afterbody 
based  on  crossflow  Mach  number,  and  an  axial  position  of 
separation  was  chosen  for  the  body  alone  so  that  the  normal 
force  was  correctly  predicted.  The  pitching  moments  of  the 
body  alone  and  body-tail  were  then  predicted  from  the  theory 
using  the  theoretical  center-of-pressure  positions.  The 
resulting  theoretical  characteristics  of  the  store  SMF 
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are  shown  in  Figure  37  for  =  0.6  together  with  the 
measured  normal- force ,  pitching-moment,  and  rolling-moment 
characteristics.  The  rolling  moment  is  associated  with 
2°  of  cant  of  all  four  tail  fins.  It  is  seen  that  the 
use  of  the  experimental  data  in  this  manner  yields  good 
theoretical  estimates  of  the  body  and  body-tail 
characteristics.  The  computer  program  should  therefore  not 
give  spurious  store  loads  for  the  various  configurations 
because  of  inaccuracies  in  predicting  store  alone 
characteristics  in  a  parallel  flow. 

In  Figure  38,  the  loads  on  store  SMF  in  combination 
with  the  clean  airplane  at  Zp/D  =  0  are  shown  for  the 
angle  of  attack  range  from  0°  to  16°  and  compared 
with  the  predictions  of  the  computer  program.  The  normal 
force  is  predicted  well  up  to  about  8°,  and  is  about 
half  of  its  free-stream  value  for  the  same  angle  of  attack. 
The  side  force  is  predicted  well  up  to  about  12°.  The 
pitching  moments  and  yawing  moments  are  not  well  predicted, 
but  they  are  generally  less  than  those  of  Figure  37  for  the 
store  in  the  free  stream.  The  center-of-pressure  position 
for  yawing  moment  is  in  error  by  as  much  as  0.5  diameters. 

In  the  theory,  the  separation  position  on  the  boattail 
corresponding  to  ag  has  been  used.  The  results  of  Figure  36 
show  that  the  average  angle  of  attack  on  the  boattail  is 
much  less  than  .  Accordingly  separation  occurs  more  aft 
than  predicted.  This  yields  more  download  on  the  boattail 
and  more  nose-up  moment.  However,  the  tail  will  gain  upload 
as  a  result  of  decreased  separation  resulting  in  a  nose-down 
moment.  What  we  observe  is  the  net  effect  of  those  two 
opposing  tendencies.  It  thus  appears  that  for  engineering 
purposes  an  improvement  to  the  computer  program  can  be  made 
for  boattail  bodies  by  taking  account  of  the  local  angle  of 
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attack  and  sidewash  on  the  boattail  in  specifying  the 
separation  location. 

The  rolling-moment  coefficient  is  small,  and  is  fairly 
well  predicted. 

Similar  results  in  Figure  39  for  Zp/'b  =  1.0  exhibit 
similar  results  as  for  Zp/D  =  0. 

COMPARISON  BETWEEN  DATA  AND  THEORY 
FOR  EFFECT  OF  THE  PYLON 


The  effects  on  the  flow  field  of  adding  the  pylon  to  the 
airplane  can  be  determined  by  subtracting  data  for  configura¬ 
tion  1  from  that  for  configuration  2.  In  this  way  the 
effects  of  the  pylon  can  be  examined  unmasked  by  configuration 
1  effects.  The  next  four  figures  show  the  pylon  effects  on 
upwash  and  sidewash  at  =  0.6  for  the  four  following 
conditions : 


Figure  as  Zp/D 


40  0  0 

41  0  1 

42  4°  0 

43  4°  1 


It  is  clear  that  the  pylon  effect  on  the  flow  field  is 
small  and  accurately  predicted  for  the  conditions  shown. 


It  would  be  normally  assumed  that  if  the  flow  angles 
induced  by  the  pylon  at  the  store  location  are  small,  that 

172 


<  / 


(c)  Pitching  moment. 


FIGURE  39.  Concluded. 


NWC  TP  6210 


- r 

- 1 - 1 - r 

_ _ _ n _ 

1 

^  o 

1 

U  ^ —  THEORY 

_ 1  1  1  1 

3  -1 

-2  -3  -4 

PROBE  AXIAL  POSITION,  X 

(a)  Upwash  angle. 

-5  -( 

P,  IN. 

-THEORY 


PROBE  AXIAL  POSITION,  Xp,  IN. 
(b)  Sidewash  angle. 


FIGURE  40.  Effect  of  Adding  a  Pylon  to  the  Clean 
Airplane  on  the  Flow  Angles  Along  the  Centerline 
Position  of  the  Bottom  Store  at  M  =0.6: 
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the  corresponding  loads  induced  by  the  pylon  on  the  store 
would  be  small.  As  an  example,  consider  the  following  table 
which  shows  the  loads  on  the  SMF  store  for  the  four  above 
cases  due  to  the  addition  of  the  pylon. 

TABLE  8.  Loads  Due  to  Pylon. 


(a)  as  =  0°  Zp/D  =  0 


Zp/D 

ACN 

ACY 

ACLM  ACLN 

ACLL 

Data 

0.07 

0.063 

-.019 

-.206  0.028 

0.002 

Theory 

0 

.031 

-.010 

-.145  .001 

0 

(b)  ag  =  0°  Z 

p/D  =  1.0 

Data 

1.20 

0.030 

-.008 

-.100  0.012 

0.001 

Theory 

1.00 

.  016 

-.004 

-.077  .001 

0 

(c)  a  =  4°  Z 

s 

p/D  =  0 

Data 

0.10 

0.069 

-.061 

-.220  0.071 

0 

Theory 

0 

.  028 

-.048 

-.138  .051 

0.001 

(d)  ag  =  4°  Z 

p/D  =  1.0 

Data 

1.24 

0.031 

-.024 

-.096  0.031 

-.003 

Theory 

1.00 

.020 

-.022 

-.072  .029 

0 

Generally  speaking  the  difference  between  experiment  and 
theory  is  not  la’ge  so  that  the  loads  were  satisfactorily 
pred ic  teu . 

An  interpretation  of  the  above  loads  will  now  be 
attempted.  The  largest  change  in  CN  is  0.069  and  the 
largest  change  in  CLM  is  -.220.  Based  on  the  slopes  of  the 
CN  and  CM  curves  at  the  origin  in  Figure  37,  those  values 
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correspond  to  free-stream  angles  of  attack  of  the  store  of 
1°  and  2.2°,  respectively.  No  such  large  pylon-induced  flow 
angles  are  to  be  found  in  Figures  40-43.  The  combined 
effects  of  flow  curvature  and  mutual  interference  appear  to 
have  the  possibility  of  magnifying  the  effects  over  what 
might  be  expected  on  the  grounds  of  induced  flow  angle  alone 
when  the  store  is  in  close  proximity  to  the  interfering 
component . 


COMPARISON  OF  DATA  AND  THEORY  FOR 
THE  EFFECT  OF  ADDING  THE  TER 

The  effects  on  the  flow  field  of  adding  the  rack  to  the 
airplane  pylon  combination  are  shown  in  Figures  44-47. 

These  results  for  both  data  and  theory  represent  configuration 
3  minus  conf iguration  2.  The  first  thing  that  is  clear  is 
that  the  induced  angles  of  upwash  and  sidewash  are  generally 
several  times  larger  than  those  due  to  the  pylon  at  Zp/D  =  0. 
It  is  also  apparent  that  the  effect  of  the  rack  attenuates 
much  in  going  from  Zp/D  =  0  to  Zp/D  =  1.0.  For  instance, 
the  induced  upwash  maxima  fall  from  about  -0.04  to  about  -0.01 
in  this  distance.  It  happens  that  the  Zp/D  =  0  position  is 
about  one  diameter  below  the  centerline  of  the  rack  and  the 
Zp/D  =  1.0  position  is  about  two  diameters  below  the  rack 
centerline.  This  suggests  that  the  upwash  is  strongly  source 
dominated.  The  effect  of  the  rack  is  modeled  solely  by 
sources  in  the  cheory.  However,  the  theoretical  source 
effect  is  too  weak,  particularly  at  Zp/D  =  0. 

Looking  at  the  sidewash  angle,  the  maximum  sidewash 
exists  at  Zp/D  =  0  and  attenuates  greatly  at  Zp/D  =  1. 

The  theory  predicts  no  sidewash  since  the  rack  sources 
produce  none  directly  below  themselves.  It  is  clear  that 
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FIGURE  46.  Effect  of  Adding  TER  to  Pylon  on  the  Flow 
Angles  Along  the  Centerline  Position  of  the  Bottom 
Store  at  M  ■  0.6;  a_  =  4°:  z_/D  =  0. 
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FIGURE  47.  Effect  of  Adding  TER  to  Pylon  on  the  Flow 
Angles  Along  the  Centerline  Position  of  the  Bottom 
Store  at  M.,  =  0 . 6 ;  ag  =  4° ;  Zp/D  =  1.0. 
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mechanisms  to  produce  sidewash  due  to  addition  of  the  rack 
must  be  introduced  into  the  computer  program. 

It  might  appear  that  the  introduction  of  doublets  to 
cancel  both  upwash  and  sidewash  velocities  along  the  rack 
are  called  for  since  they  are  not  used  in  the  computer 
program.  However,  it  appears  probable  that  doublet  effects 
will  attenuate  too  fast  to  achieve  the  desired  magnitude  of 
effect.  If  one  compares  the  effects  of  the  pylon  in  Figures 
40-43  with  the  effects  of  the  rack  in  Figures  44-47,  one  is 
struck  by  the  similarity  between  the  qualitative  behaviors 
of  the  two.  The  rack  acts  like  a  pylon  not  a  body  of 
revolution.  It  thus  appears  that  the  small  “rack  pylon"  can 
be  modeled  as  an  extension  of  the  main  pylon.  The  rest  of 
the  rack  can  still  be  modeled  as  a  body  of  revolution  with 
volume  since  it  will  not  change  the  pylon  normal  velocity 
boundary  conditions. 


It  is  not  to  be  expected  that  the  loads  on  store  SMf.  due 
to  addition  of  the  rack  to  the  pylon  will  be  predicted 
well  for  Zp/D  =  0  since  the  upwash  and  sidewash  are 
significant  at  this  location,  and  the  theory  underpredicts 
them.  It  is  of  interest  to  examine  the  load  increments  due 
to  the  rack  in  a  similar  form  to  that  for  the  pylons. 


TABLE  9.  Loads  Due  to  Rack. 


c 

MF' 


M 


OO 


0.6 


Data 

Theory 


(a)  as  =  0°  Zp/D  =  0 


Zp/D 

ACN 

ACY 

ACLM 

ACLN 

ACLL 

0.07 

0.116 

-.059 

-.222 

-.003 

0.007 

0 

.051 

.  001 

-.035 

-.002 

0 
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TABLE  9.  (Contd.) 


(b>  as 

=  0°  Zp/D  =1.0 

Zp/D 

ACN 

ACY  ACLM 

A  CLN 

A  CLL 

Data 

1.21 

0.037 

-.015  -.086 

0.045 

0.001 

Theory 

1.00 

.024 

.001  -.027 

-.001 

0 

(O 

=  4°  Zp/D  =  0 

Data 

0.07 

0.042 

-.196  0 

-.034 

+  .039 

Theory 

0 

.  050 

+.001  0.031 

-.001 

0 

(d) 

=  4°  Zp/D  =  1.0 

Data 

1.23 

0.019 

-.043  -.063 

0.040 

0 

Theory 

1.00 

.024 

0  -.024 

0 

0 

Examination  of  the  foregoing  table  immediately  shows 
that  the  computer  program  gives  no  contribution  to  CY,  CLN, 
and  CLL  due  to  addition  of  the  rack.  This  shortcoming  of 
the  method  is  more  important  at  ots  =  4°  than  as  =  0° 
since  the  sidewash  is  greater  at  as  =  4°.  An  error 
of  about  0.2  in  side- force  coefficient  occurs  and 
about  0.045  in  yawing-moment  coefficient.  There  is  a 
surprising  effect  on  rolling-moment  at  =  4°  and 
Zp/D  =  0,  the  value  of  CLL  for  configuration  3  being 
0.064  and  for  configuration  2  being  0.025.  There  is  a 
sidewash  gradient  between  the  top  and  bottom  fins  of  the  X 
arrangement  which  is,  if  anything,  weaker  at  Zp/D  =  0 
than  at  Zp/D  =  1.0.  Since  the  effect  does  not  occur  at 
Zp/D  =  1.0,  some  other  phenomenon  must  be  producing  this 
difference  which  is  equivalent  to  about  2°  of  cant  of 
all  fins.  The  phenomenon  is  believed  due  to  a  trailing 
vortex  from  the  rack  pylon. 
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The  magnitudes  of  the  largest  measured  changes  in  Tables 
8  and  9  are  now  compared. 


(ACN) 

max 

(ACY) 

max 

(ACLM)max 

(ACLN)  maJ£ 

( ACLL) 

max 

Table 

8 

0.069 

0.061 

0.220 

0.071 

0.001 

Table 

9 

.116 

.  196 

.220 

.045 

.039 

The  maximum  effect  of  the  rack  on  forces  and  rolling  moment 
is  much  greater  than  that  of  the  pylon,  but  its  maximum 
effect  on  pitching  moment  is  about  the  same. 

The  lower  store  is  released  with  the  two  shoulder  stores 
in  position,  but  for  diagnostic  purposes  it  is  interesting 
to  see  how  theory  and  experiment  compare  for  the  way  the 
store  loads  vary  with  Zp/D  when  the  shoulder  stores  are  not 
present.  This  information  is  shown  in  Figure  48  for 
Mro  =  0.6  and  otg  =  4°.  The  theory  and  experiment  both  show 
wiggles  as  Zp/D  approaches  zero.  However,  in  the  last  0.1 
a  large  change  in  the  data  to  the  attached  load  is  shown. 

The  theory  also  shows  a  sharp  change  but  in  the  opposite 
direction.  The  side-force  coefficient  goes  smoothly  to 
Zp/D  =  0  in  both  theory  and  experiment,  although  the  theory 
is  inaccurate  near  the  rack  for  reasons  already  pointed  out. 

The  data  for  the  pitching-moment  shows  an  unexpected  and 
unexplained  phenomenon  in  its  os1  illation  near  Zp/D  =  0. 

It  was  thought  that  this  result  .night  be  due  to  spurious 
data  so  comparable  results  for  =  0.8  are  shown  on 
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(b)  Side  force. 


FIGURE  48.  Continued. 
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STORE  VERTICAL  POSITION,  Zp/D 
(e)  Rolling  moment. 
FIGURE  48.  Concluded. 
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Figure  48(c) .  The  phenomenon  is  shown  but  at  a  much 
reduced  magnitude.  This  phenomenon  disappears  when  the 
shoulder  stores  are  present. 

The  yawing  moment  and  rolling  moment  show  a  change  in 
their  trends  with  Zp/D  for  small  values  of  this  parameter. 

For  the  yawing  moment  this  is  readily  ascribed  to  rack- 
induced  sidewash  effects  neglected  in  the  computer  program, 
but  the  rolling-moment  results  are  not  so  easy  to  explain. 
Referring  to  Figure  5,  it  is  seen  that  the  tail  fins  are 
well  behind  the  rack  and  the  rack  pylon.  The  possibility  of 
a  trailing  vortex  from  the  rack  pylon  seems  the  only 
possibility  which  might  induce  a  rolling-moment  coefficient 
change  from  0.04  to  0.09.  The  rack  induced  effect  is 
larger  than  that  due  to  the  2°  cant  on  the  fins. 

COMPARISON  BETWEEN  EXPERIMENT  AND  THEORY  FOR 
EFFECTS  OF  ADDING  THE  SHOULDER  STORES 

The  effect  of  adding  the  shoulder  stores  will  be 
determined  by  taking  the  difference  between  the  results 
for  configuration  5  and  configuration  3.  Any  effects 
of  the  fins  of  the  shoulder  stores  will  not  be  included 
by  this  means.  Also,  the  theory  does  not  include  the 
effect  of  shoulder  store  fins,  and  these  have  been  shown  to 
b  t-'  small. 

Figures  49-52  show  the  effects  on  the  upwash  and  sidewash 
angles  along  the  store  axis  position  for  several  angles  of 
vr.irk  and  vertical  positions  at  M  =  0.6.  Both  experiment 
and  theory  are  shown.  The  upwash  angles  get  to  be  nearly 
as  large  as  -6°,  while  the  sidewash  is  generally  much  less 


193 


0.04 


0.02 


-0.02 


-0.04 


0.02 


PROBE  AXIAL  POSITION,  Xp,  IN. 
(a)  Upwash  angle. 


PROBE  AXIAL  POSITION,  Xp,  IN. 

(b)  Sidewash  angle. 
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than  1°.  The  theory  in  some  cases  predicts  the  upwash  well, 
but  more  often  underpredicts  it.  The  degree  of  agreement  is 
surprising  in  view  of  the  fact  that  only  the  volume  effect 
of  the  two  stores  are  modeled  with  no  interference  between 
them.  Thus  the  theoretical  VZ/V  increment  due  to  the 
shoulder  stores  at  a  fixed  value  of  Zp/D  does  not  vary 
with  a  .  The  theoretical  sidewash  increment,  VY/V,  is  small. 
By  comparison  of  Figures  49(a)  and  51(a)  it  is  clear  that 
some  angle-of-attack  variation  of  VZ/V  due  to  the  shoulder 
stores  is  apparent. 

At  Zp/D  =  0  some  sidewash  exists  but  at  Zp/D  =  1  it  is 
negligible  for  a  =  0°.  At  =  4°,  the  sidewash  is  stronger 
at  Zp/D  =  0,  but  again  negligible  at  Zp/D  =  1.0.  The  rapid 
decay  of  this  sidewash  suggests  that  some  dipole 
distribution  is  causing  it.  The  two  shoulder  stores  are 
subject  to  sidewash  which  are  not  equal.  Dipoles  to  cancel 
these  sidewash  boundary  conditions  will  produce  differential 
sidewash  under  the  rack.  This  represents  a  possible  source 
of  the  sidewash.  Likewise  doublets  to  cancel  the  downwash 
distribution  along  the  shoulder  stores  would  modify  the 
upwash  and  could  account  in  part  for  the  differences  between 
experiment  and  theory.  Some  account  of  mutual  interference 
between  shoulder  stores  and  rack-pylon  may  be  required  to 
get  accurate  flow  fields.  The  shoulder  stores  can  change 
the  rack-pylon  lifting  surface  boundary  condition,  the  pylon 
can  influence  the  boundary  conditions  of  the  shoulder 
stores,  and  the  shoulder  stores  can  interfere  with  each 
other.  These  interferences  can  be  evaluated  to  see  which 
are  of  sufficient  magnitude  to  influence  the  flow  field. 

The  work  of  Martin  (Reference  8)  will  be  useful  in  this 
connect  ion . 
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It  is  of  interest  to  see  the  magnitude  of  the  loads  on 
the  bottom  store  induced  by  the  shoulder  stores.  The 


following 

table  presents  the 

resuls  for  two 

angles 

of  attack 

and  two  vertical  positions. 

TABLE  10. 

Loads  Due 

to  Shoulder  Stores. 

SMF ; 

M  = 

00 

0.6 

(a)  as  = 

0° 

Zp/D  =  0 

Zp/D 

ACN 

ACY 

A  CLM 

A  CLN 

ACLL 

Data 

Theory 

0.07 

0 

0.456 
.  334 

-.018 
.  008 

-1.049 

-.605 

-.049 

-0.012 

0.006 

0 

(b)  a  = 
s 

0° 

Zp/D  =  1. 

0 

Data 

Theory 

1.2 

1.0 

0.201 

.190 

.008 

.003 

-.452 

-.380 

-.021 

-.010 

0.003 

0 

(c)  as  = 

4° 

Zp/D  =  0 

Data 

Theory 

0.07 

0 

0.237 
.  320 

.  063 
.006 

-.955 

-.551 

-.022 

-.008 

0.028 

-.005 

(d)  as  = 

4° 

Zp/D  =  1. 

0 

Data 

Theory 

1.23 

1.0 

0.130  0 

.181  0 

.  005 

-.389 

-.385 

-.070 

-.007 

0.013 

0 

The  table  shows  that  the  addition  of  shoulder  stores  lias 
its  maximum  effects  on  CN  and  CLM,  and  has  generally  small 
effects  on  CY,  CLN,  and  CLL.  This  result  is  in  general 
accordance  with  the  flow-field  comparisons.  It  is  noted 
that  the  changes  in  CN  and  CLM  are  better  predicted  at 
Zp/D  =1.0  than  at  Zp/D  =  0,  a  fact  also  in 
agreement  with  the  flow-field  results.  Some  error  at 
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Zp/D  =  0  is  due  to  prediction  of  the  flow  field.  But 
even  with  an  accurate  flow-field  prediction,  some  errors  in 
loads  will  be  predicted  at  this  position.  It  is  probable 
that  mutual  interference  between  all  three  stores,  at  least 
for  volume  effects,  will  influence  the  loads.  This  could  be 
investigated  by  the  method  of  Martin  (Reference  8). 


COMPARISON  BETWEEN  EXPERIMENT  AND  THEORY 
FOR  EFFECTS  OF  MACH  NUMBER 

In  testing  the  computer  program  for  its  ability  to 
predict  Mach  number  effects,  we  will  compare  measured  and 
predicted  changes  in  flow  angles  and  store  loads  between 
M  =0.6  and  M  =  0.95.  This  will  be  done  for  the  clean 

CO  oo 

airplane  configuration  so  that  any  problems  the  computer 
program  may  have  predicting  the  interference  among  pylon, 
stores,  and  rack  at  M^  =  0.6  will  not  cloud  the  comparisons. 

The  normal-force  and  pitching-moment  curves  of  store 
SMF  in  the  free  stream  at  Mto  =  0.95  are  input  into  the 
computer  program  for  the  load  comparisons.  The  method  for 
generating  the  theoretical  characteristics  so  that  they 
best  fit  the  experimental  data  has  already  been  described. 
The  fit  between  experimental  and  theoretical  store-alone 
characteristics  are  shown  in  Figure  53.  It  is  seen  that 
the  data  for  the  finned  store  are  matched  fairly  well  except 
for  pitching-moment  coefficient  at  a  =  16°.  Hence  the 
body-alone  pitching-moment  curve  is  underestimated  in  the 
high  range  of  angle  of  attack.  It  is  not  possible  to  match 
the  data  for  all  four  curves  precisely  with  a  fixed 
separation  point  location  even  though  this  was  well 
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T  T 


approximated  for  the  M^  =  0.6  data  in  Figure  37.  Since  we 
are  interested  in  stores  with  fins,  the  emphasis  must  be 
placed  on  fitting  the  data  for  this  case. 

Figure  54  shows  the  afterbody  separation  locations 
used  in  the  computer  program.  The  separation  location 
moves  slightly  forward  with  increase  in  angle  of 
attack.  The  effect  of  Mach  number  is  small. 

The  comparisons  between  experiment  and  theory  are  shown 
for  upwash  and  sidewash  angles  in  Figure  55  for  ag  =  0°. 

In  this  figure  the  differences  in  upwash  and  sidewash  angles 
between  =  0.95  and  =  0.6  are  shown  along  the 
centerline  position  the  store  would  occupy  if  mounted  on  the 
bottom  station  of  the  rack  on  the  left  inboard  pylon.  The 
differences  as  measured  and  as  predicted  for  upwash  are 
small  up  to  about  a  =  12°  and  at  a  =  16°  errors  of 

b  b 

the  order  of  1°  to  1.4°  exist  in  predicting  the  Mach 
number  effect  on  the  change  in  upwash  angle.  The  agreement 
with  regards  to  sidewash  angle  is  good  almost  to  8° 
where  it  is  fair,  errors  of  the  order  of  0.5°  to 
0.75°  occurring  here.  However,  in  lieu  of  data,  the 
predictions  can  probably  be  used  for  estimating  the  flow 
angle  to  =  16°  at  M^  =  0.95  for  preliminary  design 
purposes  since  the  upwash  and  sidewash  angles  with  respect 
to  the  store  are  fairly  large  under  this  condition. 

The  differences  in  the  store  loads  between  M,  =  0.95 
and  M  =  0.6  as  measured  and  as  predicted  are  given  in 
Figure  56  as  a  function  of  angle  of  attack.  These  loads 
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FIGURE  55.  Differences  Between  Flow  Angles  at 
M.  =  0.95  and  M  =0.6  for  Clean  Airplane 
'Configuration  1)  at  Centerline  Position 
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FIGURE  55.  Continued 
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(c)  Pitching-moment  coefficient. 

FIGURE  56.  Differences  in  Loads  on  Store  Sjvip 
Between  M,,-  =  0.95  and  0.6  as  Measured  and  as 
Predicted;  Configuration  1,  Zp/D  =  0. 
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correspond  to  store  SMF  at  the  Zp/D  =  0  position.  In 
general  the  loads,  with  the  exception  of  yawing  moment, 
do  not  show  much  change  due  to  Mach  number  between 
ot  =  0°  and  8°.  As  the  angle  of  attack  increases,  the 
agreement  between  prediction  and  data  generally  becomes 
worse.  It  thus  appears  that  the  compressibility  effects  on 
loads  for  ag  >  8°  at  =  0.95  are  not  well  predicted  by 
the  linear  theory.  Some  of  this  limitation  of  linear  theory 
is  due  to  inaccuracies  in  predicting  the  flow  angles  in  the 
high  ang le-of-attack  range  as  shown  in  Figure  55. 


CONCLUSIONS 

ATTACHED  VERSUS  GRID  LOADS 

1.  The  forces  and  moments  on  a  store  in  close  proximity 
to  a  TER  can  exhibit  large  changes  within  the  first  few 
tenths  of  a  store  diameter  from  the  attached  position. 

2.  It  is  not  general  1 y  Feasible  using  present  measurinq 
methods  to  extrapolate  loads  measured  on  a  CTS  supported 
model  to  those  for  the  attached  position. 

3.  Strong  aerodynamic  interference  forces  act  on  a 
store  in  close  proximity  to  a  TER  which  is  only  weakly 
dependent  on  Mach  number  in  the  range  0.6  M  0.95. 
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ATTACHED  LOADS 

1.  Adding  the  shoulder  stores  to  the  rack  had  the 
following  effects  on  the  attached  loads  of  store  S^p  at 

=  0.6. 

(a)  The  normal  forces  received  a  large  positive 
increment  due  to  addition  to  the  rack  of  the 
outboard  shoulder  store  in  the  low  angle  of 
attack  range,  and  further  increments  of 
comparable  magnitude  by  subsequent  addition  of 
inboard  shoulder  store.  At  high  angles  of  attack, 
otg  >  10°,  the  increments  became  much  smaller. 

(b)  The  changes  in  side-force  coefficient  were  small. 

(c)  There  were  definite  changes  in  the  pitching-moment 
coefficient,  but  these  were  greater  than  those  in 
yawing-moment  coefficient. 

(d)  Rolling-moment  coefficient  exhibited  significant 
changes . 

2.  Adding  fins  to  the  shoulder  stores  under  the 
conditions  of  (1)  generally  causes  small  changes  in  the 
store  loads  except  at  large  angles  of  attack. 

3.  At  M„,  =  0.95  the  same  quantitative  configuration 
effects  described  in  (1)  and  (2)  were  found. 

4.  The  effects  of  Mach  number  on  the  attached  loads  for 
store  SAp  mounted  on  configuration  6  were  small  for  the 
range  0.6  <  Mv  <  0.95  except  for  normal-force  coefficient. 
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The  effect  of  on  CN  was  systematically  to  reduce  it  by 
an  amount  which  increased  as  the  angle  of  attack  reduced 
from  positive  (16°)  to  negative  (-4°).  These  results 
suggest  that  systematic  interference  effects  were  present. 

5.  At  Mw  =  0.6  addition  of  the  fins  to  the  shoulder 
stores  causes  negligible  effect  on  the  loads  of  store 
SAF  except  for  CLM  and  CLL  at  large  angles  of  attack. 

6.  The  addition  of  the  inboard  shoulder  store  adds 
increments  to  CN,  CLM,  and  CLN  which  reverse  sign  with 
increases  in  angle  of  attack.  It  causes  no  change  in  CY  and 
decreases  CLL.  This  behavior  is  representative  of  the  = 
0.6  and  0.95  results. 

GRID  LOADS 

1 .  At  =  0.6  and  at  us  =  0°  with  store  SMF 
within  1  store  diameter  of  the  rack  the  most  important 
things  influencing  the  bottom  store  loads  are  in  decreasing 
order  of  importance: 

(a)  Addition  of  shoulder  stores  to  TER. 

(b)  Add  it  lot',  of  TER  to  pylon. 

(c)  Addition  of  fins  to  the  shoulder  stores. 

Changing  the  angle  of  attack  to  8°  does  not  change  the 
foregoing  conclusion. 

2.  The  addition  of  fins  to  the  shoulder  stores  for  the 
ranges  of  M:o  and  i  of  the  tests  has  a  generally  small  effect 
and  tie  effect  is  largest  at  small  values  of  Zp./D.  It 
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appears,  therefore,  that  the  addition  of  the  fins  may  be 
neglected  for  trajectory  purposes  but  not  for  accurate 
attached  loads. 

3.  The  effect  of  Mach  number  on  the  loads  of  store  SMF 
for  configuration  6  is  evident  for  all  the  coefficients. 
However,  its  effects  on  CN,  CY,  and  CLL  seems  to  be  smaller 
than  its  effects  on  CLM  and  CLN.  Some  qualitative  changes 
in  CLN  due  to  Mach  number  occur  at  the  highest  Mach  number 
(Mu  =  0.95)  for  small  Zp/D. 

4.  The  data  indicate  that  the  shoulder  stores  influence 
the  loads  on  the  bottom  store  for  about  5  store  diameters 
from  the  rack. 

5.  Sway  braces  have  negligible  effects  on  the  bottom 
store  loads  except  for  CN  and  CLM  in  immediate  proximity 
to  the  rack,  where  the  effects  are  small. 

FLOW  FIELD  MEASUREMENTS 

1.  Adding  the  pylon  to  the  airplane  has  little  effect 
on  the  flow  field  in  the  position  to  be  occupied  by  the 
attached  lower  store,  but  adding  the  rack  has  a  much 
larger  effect. 

2.  Adding  a  pair  of  shoulder  stores  to  the  rack  causes 
large  changes  in  upwash,  but  the  changes  in  sidewash  angle 
are  moderate. 
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3.  Addition  of  the  fins  to  the  shoulder  stores  has  a 
small  effect  on  upwash  behind • the  fins,  and  a  small  effect 
on  sidewash  in  front  of  and  behind  the  fins. 

4.  Mach  number  causes  minor  changes  in  sidewash  and 
upwash  between  =  0.6  to  0.95  in  the  low  angle  of  attack 
range.  A  maximum  change  in  sidewash  of  1°  at  =  0.95 
was  found  over  a  small  range  of  Xp . 

5.  The  effects  of  adding  one  and  then  another  shoulder 
store  were  to  increase  upwash  in  both  cases,  but  the 
effects  on  sidewash  were  not  clearcut.  The  effects  of 
adding  each  store  were  not  equal. 


COMPARISON  BETWEEN  EXPERIMENT  AND  THEORY 

1.  By  using  some  experimental  data  in  the  computer  input 
quantities  for  the  store-alone  characteristics  in  the  free 
stream,  satisfactory  agreement  between  data  and  the  computer 
program  wore  obtained  for  the  norrna  1  -  force ,  pitching-moment, 
and  rolling-moment  curves  of  the  store  with  tail  fins  in  the 
X  roll  orientation. 

2.  For  the  clean  airplane,  the  upwash  and  sidewash 

angles  near  the  store  location  at  the  inboard  pylon  are 

predicted  adequately  for  preliminary  design  purposes  to 

u  =  8°  but  not  up  to  t  -  16°. 
s  s 

3.  For  store  S^p  m  combination  with  the  clean  airplane 
near  the  attached  inboard  position,  the  normal  force,  side 
force,  and  rolling  moment  are  well  predictid  by  the 
computer  program  up  to  a  =  8°.  However,  the  pitching 
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moment  and  yawing  moment  are  not  predicted  well  probably 
as  a  result  of  movement  of  the  separation  position  on  the 
boattail  afterbody  due  to  the  nonuniform  flow  field  and 
its  resulting  effects  on  afterbody  and  tail  normal  forces. 

4.  Adding  the  pylon  to  the  clean  airplane  produced 
downwash  and  sidewash  changes  near  the  location  of  the 
attached  store  on  the  inboard  pylon  which  were  usually 
less  than  0.5°  and  which  were  accurately  predicted  by 
the  computer  program. 

5.  The  loads  on  store  SMF  at  Zp/D  =  0  and  1  due  to  the 
addition  of  the  pylon  were  predicted  satisfactorily.  In 
magnitude  they  correspond  in  some  cases  to  a  change  in 
angle  of  attack  of  the  store  in  the  free  stream  of  1°  to 
2°. 


6.  Adding  the  TER  to  the  airplane-pylon  combination 
caused  changes  in  downwash  angle  and  sidewash  angle  at  the 
attached  store  position.  They  are  several  times  greater 
than  those  due  to  the  addition  of  the  pylon,  and  the  theory 
generally  underpredicts  these  changes.  In  fact  the  theory 
predicts  no  sidewash  changes  since  the  rack  is  modeled  by  a 
body  of  revolutions  with  volume  only  (sources  and  sinks). 

7.  The  upwash  and  sidewash  increments  due  to  the  rack 
are  qualitatively  similar  to  those  for  adding  the  pylon. 

It  appears  that  a  "rack  pylon"  will  account  for  the 
observed  effects.  It  can  be  modeled  as  an  extension  of 
the  wing  pylon,  and  the  body-of-revolution  model  of  the 
rack  still  retained. 
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8.  Generally  the  computer  program  predicts  zero  yawing- 
moment,  side- force,  and  rolling-moment  contributions  due  to 
addition  of  the  rack  since  no  sidewash  is  predicted.  The 
normal-force  and  pitching-moment  increments  are  fairly  well 
pred icted . 

9.  The  load  increments  due  to  the  addition  of  the  rack 
can  be  several  times  greater  than  those  due  to  the  addition 
of  the  pylon. 

10.  As  store  SMF  moves  toward  the  attached  inboard 
position  from  below,  the  normal  force,  side  force,  yawing 
moment,  and  rolling  moment  as  measured  show  rapid  changes 
within  the  last  0.5  diameters  of  travel.  These  rapid 
changes  are  not  predicted  by  the  computer  program.  The 
proposed  addition  of  a  rack  pylon  model  will  change  the 
predicted  loads  in  this  interval.  Better  account  of  mutual 
interference  between  rack,  pylon,  and  store  may  also  be 
required  to  predict  accurate  attached  store  loads. 

11.  The  sudden  increase  in  rolling  moment  <.  s  Zp/D 
approaches  zero  is  thought  to  be  the  result  of  a  tip  vortex 
from  the  rack  pylon. 

12.  The  pitching  moment  exhibits  an  unexplained  non¬ 
linearity  as  Zp/D  >  0  since  it  suddenly  goes  negative  and 
then  reverses  to  its  positive  value  for  the  attached 

p>'S  Ltion. 

13.  The  upwash  distribution  due  to  adding  the  shoulder 
stores  to  the  rack  is  predicted  fairly  well  at  Zp/D  =  1.0 
but  is  underestimated  ut  Zp/D  =  0.  The  use  of  doublets 

to  cancel  the  upwash  acting  on  the  shoulder  stores  could  be 
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one  source  of  upwash  not  included  in  the  theory.  Another 
source  could  be  mutual  interference  between  shoulder  stores 
and  rack  pylon . 

14.  Theory  predicts  no  sidewash  under  the  pylon  due  to 
the  shoulder  stores,  but  a  small  amount  exists  at  Zp/D  =  0 
which  varies  with  a .  It  is  thought  that  doublets  to  cancel 
the  sidewash  on  the  shoulder  stores  might  account  for  this 
small  sidewash  which  goes  away  at  Zp/D  =  1.0. 

15.  With  regards  to  the  effect  of  the  shoulder  stores 
on  the  loads  on  the  bottom  store,  large  changes  are  evident 
in  CN  and  CM  but  generally  small  changes  in  CY,  CLN,  and 
CLL . 


16.  The  changes  in  CN  and  CLM  at  Zp/D  =  1.0  due  to 
the  shoulder  stores  are  well  predicted,  but  at  Zp/D  =  0 
they  are  not  well  predicted.  While  improving  the 
predictions  of  the  flow  fields  at  Zp/D  =  0  should 
improve  this  situation,  it  may  require  accounting  for  mutual 
interference  among  all  three  stores  and  the  rack  to  estimate 
loads  accurately  in  the  attached  position  and  its  proximity. 


RECOMMENDATION  FOR  IMPROVEMENTS  TO  COMPUTER  PROGRAM 
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program  based  on  the  theory.  The  suggestions  are  aimed 
principally  at  determining  the  loads  on  the  bottom  store  in 
its  attached  position  or  close  to  the  rack.  In  solving 
this  problem  mutual  interference  between  rack  and  shoulder 
stores  is  to  be  included.  Thus  after  the  bottom  store  is 
dropped,  the  loads  on  the  shoulder  stores  should  be  more 
accurately  predicted. 

The  determination  of  the  loads  on  the  bottom  store  under 
the  conditions  of  interest  requires  that  the  flow  field  at 
the  position  to  be  occupied  by  the  bottom  store  be  accurately 
calculated  with  the  two  shoulder  stores  present.  Certain 
improvements  in  the  model  are  needed  to  achieve  this 
objective.  The  second  requirement  appears  to  be  that  the 
full  mutual  interference  among  all  three  stores  and  the  rack 
be  taken  into  account  in  determining  the  loads  on  the  bottom 
store.  Let  us  take  up  these  two  requirements  serially. 

The  present  rack  model  does  not  give  accurate  flow  fields 
because  it  does  not  account  for  the  rack  pylon.  It  is 
suggested  that  the  rack  pylon  be  modeled  as  an  extension 
of  the  wing  pylon  passing  through  the  rack.  At  the  same 
time  that  part  of  the  rack  which  does  not  include  the  rack 
pylon,  can  still  be  modeled  by  a  body  of  revolution.  If 
these  changes  do  not  give  the  desired  accuracy,  then  rack 
doublet  distributions  to  offset  the  downwash  and  sidewash 
at  the  rack  centerline  can  be  introduced. 

The  present  modeling  of  the  shoulder  stores  by  source 
distributions  does  not  give  accurate  downwash  just  under  the 
rack.  It  is  suggested  that  mutual  interference  between 
these  source  distributions  be  taken  account  of  as  well  as 
with  that  of  the  rack.  This  three-body  problem  can  be 
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handled  by  the  method  of  Martin,  Reference  8.  The  three 
bodies  already  have  a  vertical  plane  of  symmetry  so  that  the 
pylon  boundary  conditions  are  met.  If  this  does  not  give  a 
sufficiently  accurate  flow  field,  the  next  degree  of 
approximation,  which  should  be  included,  is  to  consider  the 
doublet  distributions  for  the  attached  stores  due  to  upwash 
and  sidewash. 

The  second  requirement,  when  the  flow-field  calculation 
is  accurate  enough,  is  to  take  account  of  mutual 
interference  among  all  three  stores  and  the  rack  in 
determining  the  loads  on  the  bottom  store.  This  is  first  of 
all  a  four-body  problem  in  sources  which  can  be 
approximately  attacked  by  the  method  of  Martin,  Reference  8. 
Then  the  problem  should  be  done  again  for  the  doublet 
distributions  by  the  same  general  method. 

In  addition  to  these  suggestions,  certain  others  have 
arisen  as  being  useful  under  certain  circumstances.  It  would 
be  useful  if  the  experimental  data  for  the  store  alone  in 
the  free  stream  were  used  in  the  computer  program  to  improve 
its  accuracy.  (An  attempt  in  this  direction  has  been 
previously  described  herein.)  If  a  boattail  afterbody 
occurs  on  the  external  store  with  separated  flow,  increased 
accuracy  can  be  obtained  by  accounting  for  movement  of  the 
separation  location  based  on  the  local  angle  of  attack  at 
the  afterbody  position.  In  case  rolling  moments  are 
important,  account  should  be  taken  of  the  rack  pylon  trailing 
vortex  as  it  might  induce  rolling  moments  on  the  tail  fins  of 
the  stores  . 
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A  final  requirement  would  be  to  include  the  effect  of 
shoulder  store  tail  fins  on  the  loads  acting  on  the  bottom 
store . 
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c 

c 

c 

c 


c 

c 

c 

c 

c 


poor.BAM  Rf  TCAP  (INPUT.  OUTPUT.  T  APF  s=  ,  NPUT  .  T  apt  s  =  Ou  TPuT  .  T  ARM  0  > 
program  TO  RETRIEVE  RANDOM  pun  NUMBERS  FROM  T  mF 
captive  loads  tapf 


D I MENS  ION  H(PH>  •  InDRlPS).  IPTPIIOO)  .ITITLF  ID  .RUNNO II 00) 
DIMENSION  ICONFITI.ISTORFIH 

DATA  ICAPD  /4HC»PT/.  I T APE  /*H  CAP/ 


READ  TAPE  TITlF  AND  RUN  NOS. 


RF  AD  IS.  1000,  ITITLF 
IF  (ITITLF  III  .FO.  1CARD)  0,0  TO  ?0 
WRITE  (6»?000)  ( ITITLF  < I>  «I  =  1«3T 

STOP 

?0  CONTINUE 

RFAO  (10.1010)  IHI1R 
IF  ( I  MOP ( ? )  .FO.  ITAPF)  00  TO  AG 
WRITE  ( 6 •  ?0  1  0 )  (  THOR  (  I  )  .  1  =  1  •  ?*>> 

STOP 

*0  CONTINUE 

READ  ( S • 1 0?0 )  N 

READ  ( S. 1 0  TO )  (RUNNO ( I  1  • 1 =1  »N) 
WRITE  (6.?0?0)  ( I .RUNNO < I) • 1 = 1 »N> 


C 

C 

c 


c 

C 

c 


c 

c 

c 


c 

C 

c 


SORT  THE  RIJN  NOS.  AND  THF  POINTERS  TO  The  run  NOS. 


DO  SO  1=  1 .N 
!PTR(I>=  I 
SO  CONTINUE 

IF  (N  .LF.  1)  GO  TO  AS 
N1=  N-l 
00  AO  1  =  1  .N l 
11=  1*1 
no  80  j=  7 1 «  n 

IF  (RUNNO ( IPTR ( I ) I  .LF.  RUNNO ( I PTR ( J ) ) )  GO  TO  AO 
I TEMP  =  IPTR ( I ) 

IPTR ( I ) s  IPTR ( J) 

IPTR  (  I)  -  I  TFMP 
AO  CONTINUF 
AS  CONTINUE 

SFARCH  TAPF  FOR  SORTED  PUN  NOS. 


READ  (10.10A0)  (All) .1=1, ?G I 
OLPRUN=  0.0 
DO  ISO  ! =  l.N 

no  while  desired  Run  no.  (.pfater  than  tapf  pun  no. 

100  CONTINUF 

IF  (RUNNO  (  IPTP  I  I  I  )  .1  I.  All!)  GO  TO  IhO 

IF  ( A ( 1 )  ,NF .  RUNNO ( IPTR ( I ) ) )  GO  TO  1*0 
JF  l  A  l  I  I  •  F  0 .  Oi  DR'  IN )  GO  TO  I  /A 

HPIff  Kuni'f,  FOP  Nf  W  RUN  NO. 


f  A  (  l  rONFlU  |l.(|AI.  I  C0*|F  .  I  S  T  DPF  ) 
WP  I  T  F  IP.  3010)  ()hOp(  I).  1=1"  ?SI 


I  hi  in-  Jill 
*p  I  TF  (  ».  .  S04  0  l 
wp  I  tf  i  h,  ir.sn  i 


I  Pi  IN.  I  irONf  (  Jl.»  J= 1  .  T  )  •  (  I  STORE  I  J  !  •  J=  1 

(JIJ1  .  I  a  l  J  )  «  1  =  *■  .  J  '  .'MS)  .  G(  R ) 


1) 


?. 

1. 


F> . 

7  . 
A. 
R. 
10. 

1  1  . 
1?. 
n. 

i*. 

is. 

ife. 

IT. 

I  A. 
\9. 
?0. 
?\ . 
??. 
PI. 
?4  . 
?S. 

?7. 

?B. 

PR. 

10. 

II  . 

1?. 

13. 

14. 

15. 

IS. 

IT. 

1 A  . 
1R. 
40  . 
41. 
4?. 
41. 
44. 
4  S  . 
4f>  • 
4  T  . 
4  A  . 
4R. 

50. 

si  . 

S?. 

51. 

54. 

55. 
Sf>. 
ST. 
SA. 
SR. 
SO. 
HI  . 

>'/>l 

Hi. 
H4  . 
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*ritf  < <s . o > 
oldrun*  m  1 1 1 
120  CONT  ini  if 

write  data 

WRITE  (G.2070)  R<21 ) . (H(J) .J=l 1.20) . (P(J> .J=?2.2S> 
140  CONTINUE 

READ  (10.1040)  IR(J> .J=1.2S> 
tF  ( F  OF (10)  )  ?on.ioo 

c,o  to  ioo 

1 60  CONTINUE 

FNO  OF  FILE  OP  ALL  DESIRED  RUN  NOS.  FOUND 

200  CONTINUE 
REWIND  10 


♦  ♦  RF AO  FORMATS 

1000  FORMAT  (1A4) 

1010  FORMAT  (2SA4) 

10?0  FORMAT  (IS) 

1010  FORMAT  (RF10.0) 

1040  FORMAT  (1P10F12.S) 

c 

C  ♦  ♦  WRITE  FORMATS 

c 

2000  FORMAT  (1M0.SK.17H*  •  •  FRROR  •  •  *. 

1  4GHINPUT  CARO  OOFS  NOT  SPFCIFV  CAPTIVE  LOADS  DATA  / 

2  SX,  llMCARD  TITLFr.  1A  4  ) 

2010  FORMAT  ( 1 h  0 «  SX.17H*  *  »  FRROR  »  «  », 

1  4GHTAPE  TITLE  OOFS  NOT  SPECIFY  CAPTIVE  LOADS  DATA  / 

2  SX.11HTAPF  TITlEs*  2SA4) 

•2020  FORMAT  ( 1H1 , 1 7HTNPUT  RUN  NUMRFBS  / 

1  2X . INI , IX , 7MRUN  NO.  /  ( 1 X, 12.F 1 0.2) ) 

2030  FORMAT  (1M1.1X.2SA4) 

2040  FORMAT  ( I  HO . 3X . 7HPUN  NO . . 7X . 20HPARENT  CONFIGURATION. 

1  1 SX . SHSTORE ,  / 

2  SX.I3.0X.7A4.7X.1A4) 

20S0  FORMAT  (1H0.17X.14H*****  FRFE-STRFAM  CONDITIONS  •*«•*.  / 

1  4X.4HMACH.4X.2HPT ,gx,2hTT.7X.1HP,7X.1HT,7X.)MV.7X, 

2  1H0.7X.SHRF/FT.  / 

1  SX  .  2HN0.4X  .4HPSF  A  .  IX  .GHI'FG  (R>  .  IX  .  4 HPSF  A  .  IX  . 

4  GHRFG <B) . 2X .GHFT/SFC . 4X , 1HPSF ,GX ,GHFT (- 1 )  / 

5  2X.FG.3.6FR.1  .F7.2.GH* 1 0**G) 

20GO  FORMAT  ( 1H0.S0X ,22H*»***  MODFL  DATA  •••»•  / 

1  4X .SHALPHA . IX .SHALPHA , IX .4HPF T A  / 

2  SX  ,  1HPAR,  4X  .SHSTORF  ,  IX  .SHSTO'RF  .4  1. 2HCN.RX  . 2MCY . GX . 

1  1HC  AT  .SX.  1HCLM.SX  ,  1HCI  N.SX.  1HCLL.GX  ,  1HNCP.ISX  . 

4  1HYCP. 4X , 1HMNF . IX , SHMOOTN. lx .2HCR.SX . 2HVR  / 

5  SX.  1HDFG.SX , THDFG.SX , 1HDFG.S1X ,ShDI AMS.4X .SHD 1  AMS, 

G  OX,  7HLBM/SFC  ) 

20  7  0  format  (  1  X  ,  1  (  2«  .F*..  2)  .*.  (  2*  .FS.  1)  ,  >  (  2X  .F7 . 1)  .4F  7  .  7  | 

C 

c 

STOP 

FNO 


OS. 
GG  • 
G7. 
OR. 
GR. 
70. 
71  . 
72. 


74. 
7S  . 
7  G. 
77. 
7R. 
7R. 
RO. 
«1  . 
P?  . 
Rl. 
R4. 
RS. 
RG. 
R  7  . 
RD  * 
PR. 
RO. 
Rl  . 
R2. 
R3. 
R4  . 
RS. 
RG. 

RR. 
RR. 
100. 
101  . 
102. 
101. 
104. 
1  OS. 
1  OG. 
107. 
1  OR. 
1  OR. 
110. 
111. 
112. 

113. 

114. 

115. 
1  1G. 
117. 
1  IP. 

1  1  R. 

1  20. 
121  . 

1  22. 
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SUBROUTINE  CONFIC.  <n,  I  CONE . I  ST ORF  > 

POO. 

c 

?0  1  . 

c 

OFCTPHFR  THF  CONE  I GnR A  T I  ON  COOF  ANO  SFT  UH  AN  ALPHABETIC 

pop. 

c 

ARRAY  FOR  THE  rONF  K-IIRAT  TON  ANTI  T  hF  STORF 

?0  3  • 

c 

P  0  4  . 

c 

POS. 

01  MENS  I  ON  IC ( ] 0) , irONF ( 7 ) . IS ( 7) . I STORF (31 

?0h. 

OATA  ir  /4HF-4C.  4H.P1,,  4HT.IS.  4HPU).. 

4H?) F , , 

P07. 

1  4H.PI  .  4HT  .  4H?)U  .  4H(S3>.  4HU  . 

4HF 

POR. 

?  4H.PCL.  4HF.P0.  4HPIF  ,  4HTM, ( .  4HS?)F. 

4H. (S3. 

POR. 

3  4HIF  / 

P10. 

OATA  IS  /4HMK-P,  4H3.ME,  4H3.MU.  4H3.AF. 

4H3.AU. 

PI  1  . 

1  4HPR0R.  4  HE  / 

PIP. 

OATA  IRLANK  /4H  / 

?  1  3  . 

c 

?  1  4  * 

TE  (fl  .E 0.  610.0)  B=  Rl.O 

PIS. 

Mfl=  R 

?  1  6  . 

MY  =  MOO  I  MR .10) 

PI  7. 

MX=  MR/10 

PIP. 

c 

?  1  R  • 

c 

CASE  CONE  I GUR AT  1  ON 

??0. 

c 

??1. 

c 

CASE  X  =  0 

???. 

c 

??  3  . 

TE  (MX  .NF.  0)  00  TO  SO 

??4. 

00  40  1=1.7 

PPG. 

I  CONE (  I  )  =  IRLANK 

??6. 

*0 

CONTINUE 

??7. 

00  TO  ISO 

PPP. 

c 

PPR . 

so 

GO  TO  (60. 70.R0.R0. ] 00. 1 1 0. 120. 1 30. 140) .MX 

?  3  0  . 

c 

?  3  1  . 

c 

CASE 

P3P. 

c 

P  3  3  . 

(so 

CONTINUE 

P  34  • 

irONF ( i ) =  rr ( i ) 

?  3S  . 

no  i=  ?.? 

?36. 

ICONF! I) =  IRLANK 

P  3  7  • 

6S 

CONTINUE 

?3P. 

no  to  isn 

PPR. 

c 

P40. 

c 

r ASF  X=? 

P4  1  . 

c 

P4P. 

70 

CONTINUE 

?43. 

TFONF (1 ) =  IC  <  1  ) 

P44  . 

KONF|?)s  1C(M 

* 

?4S. 

no  7S  T  =  1.7 

?4t>. 

ICONF  ( I 1  =  I  HI 4NK 

?47. 

7S 

FONT  I N*  fF 

?4P. 

on  to  \so 

?4-R. 

r 

P-.0. 

r 

CASE  Trl 

?S1 . 

c 

PSP. 

00 

fOMT  f  Mir 

ppp. 

I  r  ow  ( i  i  =  inn 

PP  4  . 

T  rr^f  (  -  ,  J  -  , 

pss. 

f  o  •  r  {  \  *  *  \ 

pss. 

r  /■»  at  t  .  „  ,  * 

?S7. 

T  r  ON*  •  I  •  ■  1  -U  AS* 

ppp. 

at 

r  tnt  i%«.f 

PHR. 

TO  |  ' 

PM'. 

r 

p<  1 . 

r 

r  r  i  -  a 

?pp. 

c 

Pp  3. 

0,1 

c^nt is 

/Pa  • 

"  <-  I  •  »  .  » 

PPG. 
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iconf  i  n  =  inn  ?66. 

OS  CONTINUF  ?S7. 

Iconf <*>=  mm 

no  os  s,7  ?6s. 

IC0NF<1>=  IWLANO  ?70. 

06  CONTINUE  ?7  1  . 

oo  to  isn  ?7?. 

C  ?73. 

c  C  ASF  X  =S  ?7*. 

C  ?  7S  . 

ion  CONTINUE  ?76. 

no  10S  1=  1.4  ?  7  7  • 

I CONE  (  1  )  =  IC(II  ?7B. 

10S  CONTINUF  ?7<j. 

iconfis)=  mo)  ?«o. 

ICONF (6) =  in  10)  ?«  1  . 

IC0NF(7)=  IHLANK  ?B?. 

Oo  TO  ISO  ?«3. 

C  ?F>*. 

c  0  ASF  X  =  6  ?«s. 

C  ?«6. 

110  CONTINUE  ?«7. 

no  ns  1=  1.3  ?««. 

ICONF ( I ) ric ( I )  ?H9. 

11S  CONTINUE  ?on. 

ICONF (A) =  1C(S)  ?0] . 

ICONF <S)=  IC<0)  ? o?. 

IC0NF(6>=  IC(ll)  ?9  3. 

TCONF(7)x  IHLANK  ?9*. 

00  TO  ISO  ?OS. 

?06. 

C ASF  X  =  7  ?07. 

?OP. 

1?0  CONTINUE  ?»<). 

ICONF ( 1 ) =  IC  (1  )  300. 

I  CONE  (  ? )  =  ini?)  301. 

TCONF(3)=  in?)  30?. 

ICONFU)=  IC(  3)  303. 

I CONF  ( S )  =  I  C  (  s  )  3n*. 

1C0NF(6)=  I C ( O )  30S. 

ICONFI7):  I C ( 1 3 )  306. 

00  TO  ISO  307. 

C  30ft. 

C  CASS  *  =  H  300. 

C  3  1  n  . 

1 30  CONTINUE  311. 

no  ITS  1=  1.3  31?. 

ICONF!  n=  inn  313. 

ITS  CONTINUF  31*. 

ICONF  (  *  )  =  in  I  4  )  3  1  S . 

OO  1 36  1=  S. 7  316. 

ICONF ( n  =  m  ANK  317. 

136  CONTINUF  310. 

no  to  isn  3 1 <j. 

c  nn. 

C  c  aft  < - Q  3^1. 

c  3??. 

10(1  CONTINUF  n  1. 

ICONF  (  1  )  -  TC  (  i  )  3’*.  . 

T  CONF  (?)  *  T  C  I  -■ I  3^S . 

T  r  0 v r  (  3  )  r  I  r  I  1  u  ,  3.-6  . 
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tfonf  (u)  -  iniM 

3/7 

Tfonfiaiz  i r  t  i  7 ) 

3?8 

T  FONF t  A)  =  ] 0 ( 1  A ) 

3?9 

TFONFITIr  IRLAN* 

3  *0 

c 

3  3  1 

c 

FNO  C ASF  CONFIGURATION 

3  V 

iso 

FONT  INIIE 

313 

c 

33« 

c 

3  IS 

c 

CASE  STORE 

33* 

c 

3  37 

ISTORF ( 1 ) =  IS(1> 

338 

TSTORF (  7)  =  IRLANK 

339 

GO  TO  ( 1 70, 1  AO. 1 R0. ?00.?1 O.??O.?A0)  .  MY 

14  0 

c 

34  1 

r 

C ASF  Y  =  1 

34? 

c 

343 

1  70 

CONTINUE 

344 

T  ATORF i ?) =  IS<?) 

34  8 

GO  TO  ?A0 

34* 

c 

347 

c 

C  ASF  Y  =  ? 

348 

r 

14  9 

1  AO 

FONT  TNdF 

ISO 

1ST  ORF ( ? ) =  IG(A) 

381 

GO  TO  put) 

38? 

c 

381 

c 

CASE  Y  =  3 

384 

c 

388 

1  AO 

FONT  TNUE 

186 

TAT  ORE {?) ~  I S 1  A ) 

187 

GO  TO  ?40 

388, 

c 

38Q , 

C 

C  ASF  V-U 

3*0 

c 

3*1 

?00 

FONT INUE 

3*? 

I  STORE (Pis  TS(A) 

3*3 

GO  TO  ?40 

3*4 

r 

3*8. 

c 

F  ASF  Y  =S 

3**. 

c 

3*7 

?10 

FONT  I  Ml  IF 

3*8 

ISTORF  (  1  )  =  IS!*-) 

3*9 

TSTOPE(?)=  T  A ( 7 ) 

170. 

GO  TO  ?40 

371 

r 

17?, 

c 

FAAE  Y=  A 

173, 

??o 

FONT  TNUF 

1  74  , 

T ATORF  ( ?\ -  1  A ( ?\ 

378 

GO  TO  ?U 0 

37*, 

r 

177  , 

c 

F  A  A  E  Y  =  7 

178, 

?70 

FONT INOF 

1  79, 

TAT  ORF  (  r>l  =  T  A  1  7) 

1MO  , 

r 

*8  1  . 

r 

ENO  FAAf 

W. 

c 

1H  1, 

Pun 

FONT  T  ’ (i  if 

3*4  , 

c 

1  ->8  , 

LJt  T  ,  IW\ 

i  , 

r 

3  -  7  , 

>  ✓ 
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PROGRAM  RE  T  Hun  (  INPUT  .OUTPUT  .  TAPE  1  NPUT  •  T  APE  MIIUIFU  1  .  T  APE  |  n  I 
c  PROGRAM  TO  RFERIFVF  RANDOM  PUN  NL/mHFPS  FROM  T HE  ? 

C  GPfn  TAPE 

C  * 

C  S 

DIMENSION  h  (?*)  .  |HOP  (?S)  .  IPTP(  )  On)  .  I  T  I  TLF  (1)  .PUNNOI  1  0(1)  G 

DIMENSION  ICONF ( T) . TSTORF (11  7 

C  « 

DATA  T CAPO  /4HGP in/ .  I T APF  /4H  GP 1 /  P 

C  1  0 

C  1  1 

c  PF  AO  T  APF  T I TL  F  ANO  PUN  NOS.  \? 

c  n 

PF  AO  IS.  1000)  JT1TLF  14 

IF  ( I T 1 TLF  C 1 )  .FO.  ICAR0)  GO  TO  ?0  IS 

WPITF  (G.?000)  ( ITITLF < I ) . 1=1 . 1)  1G 

STOn  17 

?0  CONT1NIIF  1« 

PE  AD  (10.1010)  IHOR  IP 

IF  ( I  HOP ( ? )  .FQ.  1TAPF)  GO  TO  0  ?0 

WPITF  (G.POIO)  ( I MDR <  1  )  •  I  =  1 . ?S  >  ?\ 

STOP  ?? 

40  FONT  I NUF  ?3 

PFAO  (S.10?0)  N  ?u 

PF  AO  (S. 10101  (PUNNOI I  I, 1  =  1 .N)  ?S 

WPITF  (G.?0?0>  < I .PUNNOI I > . 1=1 ,N)  PG 

C  ?7 

C  SORT  ThF  RUN  NOS.  AMO  T HP  POINTERS  TO  ThF  RUN  NOS.  ?H 

C 

no  so  i=  l . n  to 

IPTR(I)=  i  1 1 

GO  OONTINIIF  TP 

IF  (N  .LF.  1)  GO  TO  RS  11 

N 1 =  N-l  14 

no  ro  i=i.ni  is 

11=1*1  is 

no  80  J=  1 1 .N  17 

IF  (PUNNOI JPTB ( I ) )  .LF.  WUNNO(IPTPIJ) ) )  GO  TO  RO  TR 

ITFMP=  IPTR(I)  TP 

I P  TR (  I  )  =  IPTP(J)  40 

TPTP(,l)r  1TEMP  «1 

RO  CONTINUE  4  P 

RS  CONTINUE  4  J 

44 

SEARCH  T  APF  FOP  SORTFO  RUN  NOS.  4S 

4  G 

RF  AO  (1  0,  1  040!  (P(T)«t  =  1*  PG )  47 

OLODIIM-  o.o  4  R 

no  1  GO  I  =  1  .fi  4  P 

C  s  0 

C  OP  WHIIF  OF  S I  RF  o  RUN  NO.  GREATER  THAN  TAPE  P((N  Nil.  SI 

C  c  P 

1 00  FONT INUE  ST 

T r  (PUNNOI  TPTWI  T  )  )  .IT.  0(1))  ,1  TO  (GO  S* 

IE  (0(1)  .Ml.  PUNNOI  JPTR(  T  ))  )  (-0  TO  )40  SS 

TE  (  R  (  1  I  .10.  0|  ITPUN)  l>()  TO  )  ->0  SG 

C  S  7 

C  RRirr  hfap(ng  fop  new  pun  no.  ro 

C  so 

CALL  COllTi.  (0(10).  1  CURE  .  ISTi'Pt  I  >  (I 

W[-  I  T  E  (►.Poll))  (  I  MPH  (  I)  .  ,1  -  1  .  PS)  R  1 

|  hi  in=  •«  (  |  I  R  ■> 

WP  T  T  f  (g.POaOI  1  P|IN.  (  ]  CONE  (  J1  ,  J=  1  .  7  )  .  (  I  R  TPRE  (  II  .  .)=  1  .  T)  R( 

.  '  '  I  ‘  (  <  .  . .  I  I  (  I  >  •  I  :  P  •  .  I  .  (  P  (  I  I  ■  I  -  '  .  o  I  .  "  (  R  1  .  ■>  (  I  »  . 
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wR  i  r  >  i  <• .  ,'Hmi  t 

oi  rtwnn-  (•  i  l  t 

IPO  CONT1NIIF 

C 

c  WW I T  h  OATA 

c 

txRITF  (»i.?n?ni  M  (  11  )  .M  (  1  7»  «n  (181  >H  (  IS)  .R  (  )  61  . 
1  <R(J>.J=l?.l4).(H<Jl.J=?),Ps),R<lR).H(P0) 

140  FONTTNUF 

RF  40  (10.10401  ( H ( J>  ,  J=  1 • ?b ) 

IF  (FOF ( 1 0) I  POO. 1 00 
100  FONTINUF 

FNO  OF  FRF  OR  AU.  OF  S  1  RF  O  RUN  NOS.  FOUNO 

POO  CONTINUF 
RF  W  T  NO  10 


.  .  RF  AO  FORMATS 

1000  FORMAT  (  0  A  4  ) 

1010  FORMAT  (?SA4l 
]0?O  FORMAT  (IS) 

1000  FORMAT  ( OF  10.0) 

1040  FORMAT  ( 1R10F 1 ?.S1 

C 

C  ♦  ♦  WRITF  FORMATS 

C 

?000  FORMAT  (1H0.SX.17H4  •  •  FRROR  •  •  •. 

1  07HTNPUT  FARO  OOFS  NOT  SPF f I F  Y  GRID  DATA  / 

?  SX,  11HCAPO  TITIF:,  XA4 1 
?010  FORMAT  (  1  MO .  SX.17H*  «  «  FRROR  •  « 

1  07HTAPF  T I T  l  F  OOFS  NOT  SPF  TIFT  GRID  OATA  / 

?  SX.11HTAPF  T[TlF=.  PSA4) 

?0?O  FORMAT  ( 1H1 , 1 7HINRUT  RUN  NUMRFRS  / 

1  ?X. 1MI. T« . 7HRUN  MO.  /  (  tX.lP.FlO.PU 
POXO  FORMAT  ( 1 M 1 . 1 » . ?SA4 ! 

P040  FORMAT  (  1H0  .  OX  •  THRUM  NO . . 7 X , ?0 HP ARF NT  CONFIGURATION. 

1  1SX .SMSTORf ,  / 

?  SX.I0.8X.7A4.  7*.  (  A  4 ) 

?0S0  FORMAT  . .  FRff-STRFAM  fONOITlONS  44444,  / 

1  4X.4HMAFH.4X. ?“PT ,SX.PHTT,7X.1hP.7X,1HT.7X.1HV.7X, 

?  1H0.7X.SHRF/FT,  / 

X  SX  ,  ?HNO,  4X  .  4HHSF  A  .  IX  .  «,HIH  ,  XX  •  4  HPSF  A.  IX  . 

4  SHOFG (Rl  . ?X .SHF T/S(  F .4 X . XHPSf  .SX .SHFT  ( -1  I  / 

s  ?X  .FS.  0.SF8, 1  ,F  7.P.SH4  1  044f,| 

?0S0  FORMAT  (  1  HO. SO  X  .  ?  JH.4.44  Mi>|lFl  OA  T  A  444  44  / 

1  4X  .SHAl  PHA  .  X«  .SHAI  PHA.  XX  ,4Hff  >  A  .4  X  .4H0IHA  ,  4X  ,4HI'PHl  / 

P  SX  .  XHPAR  .  4X  .ShsTORI  ,  XX  .  mis  1  OR»  ,  (X  .SHSTORF  ,  XX  .  SMSTOHf  . 

X  4X,?HXP.SX,7HVP.4X.4H/P/0.4X,7MrM.SX.PHrY,SX, XHfAT, 

4  SX  .  XHFI.  M.sx  .  >Hfl  X  ,SX  .  u..|  l  .sx  .  XMU(  ,  OHy  FM  / 

5  SX  .  XHOFF,  ,  SX  .  4H(it  F.  ,S  x  .  XmOF  I,.  SX  ,  (••(!(  I..SX  ,  XHOF  (■’,  s  X  . 

S  PHTN.SX . ?»1M.M.1 ,S.mi AMS. 4 x .SMC  I  AMS  I 

PA 70  format  (  1  x  ,s  I  ->X  .Fs.  ;>|  .  1  (/,  ,i  s.Pl  .x-  ( /X  .Fs.  (I  ,P (  /»  .f  7.  01  1 
F 
F 

STOP 
F  Nfl 


SN. 

^1^1  * 
s  7  , 

SH  . 

SR. 
70. 
71  . 
7?. 
7  0. 


7S. 

77. 

78. 
7R. 
80. 
Ml. 
«?. 
«0. 
H4. 
MS. 
MS. 
M7  . 
MR. 
8R. 
RO. 
R 1  • 
Rp  . 
RO. 
R4. 
RS. 
RS. 
R  7  . 

RR. 
1  00. 
10). 
10P. 
1  0  X. 

1  04. 

10S. 

1  os. 
107. 
1  08. 
1  OR. 
110. 
111. 
1  IP. 
110. 
114. 
1  IS. 
US. 
117. 

I  1M. 

1  1  R. 
IPO. 
1P1 . 

1  pp. 

I  p  (. 
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Siiooo.it  Inf  roof  | r.  (<-,  ioonf  .  islnpt  i  /no. 

R01  . 

ntriPHfu  t  xi-  ro'iF  iT'Iimat  t'Vj  root  ami  st  t  up  ao  ai  o*-*AMt  t  I  r  y  ir. 

APPAY  n,o  IhF  r'lftlAiiiMI  ION  Aoo  T*<t  ST0PF 

R04  . 

/04. 

o[mfnsi<'n  i  r  (  |  >ti  .  iroNt  (  7) .  i  s  (  7 1  .  |  oToof  1 1)  /m,, 

oata  ir  /4.if  -«r .  ah. pi,.  4hT.(S.  a*M)ii..  4H7IF..  rot. 

|  AM, PI  .  4>ll  .  4HR|ll  ,  4HIS1I  .  A  HO  .  4HF  .  ROR. 

R  4H . °FL •  4HF.P0.  4HR1F  .  4HIM.I,  4HSR1F.  4H.IST.  ?A4. 

1  4H)F  .  4HN0NF  /  RIO. 

OATA  IS  /4HMK-R.  4HT.MF.  4H7.Mll.  4H7.AF.  4H7.AU.  ?||. 

1  4HPR0R,  4HE  /  R]R. 

DATA  IRLANK  /  4H  /  ?n. 

?14. 

IF  (0  .FO.  MO. 01  fl=  91. 0  RIS. 

Mp=  R  RIS. 

MY  =  MOO ( MR ,10)  R  1  7  . 

MX  =  MR/10  R1H. 

RIM. 

C  ASF  FONF  1  fil/R  AT  I  ON  RRO. 

RR1  . 

C ASF  X  =  0  RRR. 

RR7. 

IF  (MX  .NF.  01  on  TO  SO  RR4. 

irONF(l)=  IC(IO)  rrs. 

00  40  1  =  R,  7  ??T>. 

irONF(II=  I  ML  AfsIK  RR7. 

40  FONTINUF  RRR. 

00  TO  ISO  RRR. 

RIO. 

SO  00  TO  (<S0. 7O.R0.RO,  1  00.  11  0.  1/0.  1  00. 140)  .M*  R71. 

R7R. 

C  ASF  Xrl  R77. 

R14. 

00  CONTINUE  RIS. 

iroNFim  rrin  rim 

OO  SS  1=  R.7  ? 17 . 

irONFdls  IHLANK  R7H. 

6  S  rONTINIIF  ?TR. 

00  TO  ISO  R40. 

?4l  . 

C  ASF  X=R  R4R. 

?<M. 

70  rONTINIIF  R44. 

irONF(i)=  I r ( i )  / 4  s . 

TroNF(?)=  ifim  r as. 

no  7S  1=  1.7  74). 

ICONF (11=  [Ml  ANK  RaR. 

7S  rONTINIIF  /4<i. 

00  TO  ISO  /SO. 

/‘.I  . 

r ASF  x  =  7  /SR. 

/SI. 

RO  roNT  [ni  it  RM. 

I r  ONF  (  1  )  =  ir  (  1  1  RSS. 

Ifdnf  (  p  )  =  I  r  (  ?)  rso. 

IrONF  ID:  Kin  R‘-  7. 

no  QS  I s 4 . 7  RSR . 

irONF(I)=  I>'l.  AN"  ?sr. 

04  rONTJNMF  RA-0. 

OO  TO  I  4  0  R»  1  . 

r  r*-r. 
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FASF  «  =  4 


OS 


Qt)  CONTINUE 

On  os  T  =  1 . 3 

iroNFtn=  if  in 

FONT  IMIF 
T CONE (41-  I  F  <  a  ) 

HO  OS  7=  S.7 

trONF(!)=  F«LANK 
FONT [NUF 
GO  TO  ISO 


OS 


CASF  X  =  S 


100  CONTINUF 

no  70S  7  3  7.4 
irONF(7>=  I C ( I ) 
70S  FONT  7  NUF 

IF0NF<S>=  IC(O) 

7  FONF ( G I  =  7  F ( 1 0 ) 
TFONF (77=  IOLANK 
GO  TO  ISO 


F  ASF  X=G 


1  IS 


770  CONT7NUF 

no  1  IS  1=  1.3 
I FONF ( I ) = I F ( T  > 
FONT  7  NUF 
T  FONF ( 4 )  =  I F ( S ) 

7  FONF ( S 1 =  1C(01 
IFONF(S):  r  F  <  l  7  > 
T  FONF ( 7 ) =  701  AN* 
GO  TO  ISO 


F  ASF  X=7 


1?0  FONT INUf 

I FONF (17=  7  F ( 1 ) 
TFONF (?)=  I F ( 1 ? ) 
7FONF(3l=  1 F  <  ? ) 

7  FONF ( 4 )  r  1F(3) 

I FONF ( S ) =  I F  <  S I 
I FONF ( G ) =  7  F ( O ) 
TFONF (71=  7 C ( 7  3 ) 
GO  TO  ISO 


CASF  X  =fl 


1  3S 


130  F0NT7NUF 

00  1 3S  1=  1,3 
TFONF (I)=  !C(7) 
FONT INUF 
I  FONF  (  4  )  =  IFM4I 

no  i 3s  I-  s,7 

7 F ONF  (  I  )  =  JULAN"- 
FONT INUF 
GO  TO  ISO 


13* 


F ASF  »=g 


140  CONTINUF 

1 OONF  (  1  »  =  IF  I  1  I 
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!roNF(?)=  rri?) 
jrr>NFM)=  i  r  <  i  *->  > 

imMF  (4)  r  |  C  (IF) 

Iconf  is)  =  ir  ( 1 7) 

ICONF  cs>  =  jr  (  |  «i 
ICONF  <*,)•=  |r  (  |  hi 
ICONF (71=  IRLANK 

f  NO  C ASF  CONFIGURATION 

iso  continuf 


C ASF  STORF 

ISTopf (i)=  is<  1 1 
ISTOPF<3)=  IRLANK 

GO  TO  <1T0.1«O.IPO.?00.?IO.??0.?30) .  MV 

C  ASF  Y  =  1 

170  CONTINUF 

I STOPF (?)  -  I S ( ? ) 

GO  TO  ?40 

C  ASF  Y  =  ? 

l«0  CONTINUF 

istorf(?)=  ishi 

GO  TO  ?40 

C ASF  V  =  3 

ISO  CONTINUF 

I STOPF  <?) ~  IS(4) 

GO  TO  ?4  0 

C  ASF  Y  =4 

?00  CONTINUF 

I STORF  I  ? )  s  ISIS) 

GO  TO  ?4  0 

CASF  Y  =  S 

?10  CONTINUF 

I STORF (  1  )  =  ISIS) 

ISTOPF(?)=  I S ( 7 ) 

GO  TO  ?40 

CASF  Y:  S 
??0  CONTINUF 

1ST  ORE  (?)  -  IS  (?) 

GO  TO  ?40 

CASF  Y  =  7 
<*30  CONTINUF 

ISTORF(?)=  I S I  3 ) 

FNO  CASF  STORF 

7*0  CONTINUF 
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PROGRAM  RFTFI  O  (  INPUT.  OUTPUT.  T APF s= INPUT . T APF S=OUTPUT . T APF 1 0 > 

PRnr,u*M  TO  RFTRIFVF  PANOOM  PUN  NUMF>FRS  FROM  THF  ?. 

FLOW  FIFin  SMPVFY  tapf  3. 

4  . 

s. 

0  !  MENS  I  ON  R (PS) . IHORtPS) . IPTR ( 100) . ITITLf ( 3) .RUNNO! 1001  G. 

DIMENSION  ICONF (7) . TSTOPF (3)  7. 

R. 

DAT  A  1  CARD  /4HFL0W/.  ITAPF  /4H  FLO/  9. 

10. 
1  1  . 

PE  AO  TAPE  TtTLF  AND  RUN  NOS.  1?. 

13. 

READ  (S.1000)  ITITLF  14. 

IF  (ITITLFU)  .FO.  ICARD)  GO  TO  ?0  IS. 

WRITE  (6.POOO)  (ITITLF(I) *1=1.3)  10. 

STOP  17. 

?0  CONTINUE  IS. 

RE  AO  (10.1010)  THOR  IP. 

IF  (IHOR(P)  .FO.  ITAPF)  GO  TO  40  ?0. 

WRITE  (G.POIO)  ( 1 HDR (I)«I=1«?S)  PI. 

STOP  PP. 

40  CONTINUF  ?3. 

RF  AO  (S.inpO)  N  ?4. 

RF AO  (S.  1030)  (PIINNO  (  T  )  ,  1  =  1  ,N)  PS. 

WRTTF  (6.P0P0)  ( I ,RUNNO( I ) . 1=1 .N)  PG. 

P7  . 

SORT  THF  RUN  NOS.  AMO  THF  POINTERS  TO  THF  RUN  NOS.  PR. 

PR. 

00  40  1=  l.N  ‘  30. 

T  RTR ( I ) =  I  31. 

(SO  CONTINUF  3p . 

IF  (N  .LF.  II  GO  TO  FIS  33. 

N 1  =  N-l  34. 

no  RO  1=1 . N 1  3S . 

11=  1*1  36. 

00  RO  J=  1 1 .N  37. 

IF  (RUNNO ( 1PTU ( 1 ) )  .LF.  RUNMO ( I PTR ( J ) ) I  GO  TO  RO  3P. 

I TFMP-  IPTP(l)  3Q . 

IPTP ( I ) =  TPTR(j)  40. 

I PTR ( J ) =  1TFMP  41. 

RO  CONTINUF  4  ? . 

RS  CONTINUF  43. 

44  . 

SFAPCH  TAPF  FOP  SORTFO  RUN  NOS.  4S. 

44. 

RF  AH  (1  0.10401  <R (  I  1  .  I  =  1 ,PS)  47. 

OL  ORt  IN=  0.0  4  R  . 

00  ISO  [=  l.N  4  R . 

C  u  0  . 

C  00  WHlLF  OF  S 1  RF n  R(|N  NO.  GRFATFR  THAN  TAPF  RUN  NU.  SI. 

C  R?. 

100  CONTINUF  ST. 

I F  ( RUNNO  (  I P  TR  (  I  )  )  .LT.  S(II)  go  TO  Iso  S4. 

IF  (0(1)  .  N(  .  RllMNO  (  IPTR  (  I  )  )  )  GO  TO  1  4  0  ss. 

IF  (0(1)  .FO.  01  ORIIN)  GO  TO  IPO  SS. 

C  s  i. 

C  WPITF  HF  AO  I  MG  FOP  NF w  RUN  NO.  sH. 

C  SQ. 

CA|  1  CONFIG  (U(10).  ICONF  . ISTORF )  SO. 

WPITF  (G.P0301  ( I mOh ( j ) , J= 1 . ps )  si. 

1  Pi  IN=  -till  SP. 

WPITF  (S.P040I  IRON, 1 ICONF ( j) ,J=1 , 7 )  S3. 

W-ITF  (*..>0  =  01  (  »  i  I)  .  I  =  P.  „  |  .  (P  (  )l  .  J  =  S  .  H  )  (S)  .4  {  41  14, 
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wriu  (»-./>o«n>  hiiI)*h(?m  6S. 

wriif  h.posoi  66. 

OLORUNs  M ( 1 >  67. 

1P0  CONTINUE  6H. 

C  69. 

C  WRITE  0414  70. 

C  71  . 

WRITE  (6.?07n>  (MIJ) * J= 1 ?♦ 1  * ) *0(16) .  H  1  1  S  ) • R ( PS ) .  7?. 

1  (R(J) .J=17.?0) ,R<?1) .HI?*) .Of??) .R(?l)  71. 

1*0  CONTINUE  7*. 

READ  (10.10*0)  (R(J)  ..1=]  t?6> 

IF  (FOF(lOl)  POO. 100 

160  CONTINUE  77. 

7B. 

ENO  OF  EILF  OR  All.  OESIRFO  RUN  NOR.  FOUND  79. 

«0. 

POO  CONTINUE  Rl. 

REWIND  10  PP. 

P3. 

R*. 

♦  *  REAP  FORMATS  PS. 

P6. 

1000  FORMAT  (3A4)  P7. 

1010  FORMAT  (PSA*)  PR. 

10P0  FORMAT  (IS)  PR. 

1010  FORMAT  ( RE  10.0)  90. 

10*0  FORMAT  (1R10E1P.S)  91. 

C  9?. 

C  ♦  ♦  WRITE  FORMATS  91. 

C  9*. 

P000  FORMAT  (1H0.SX.17H*  •  •  ERROR  •  «  *.  9S . 

1  SOHINPUT  CARO  OOFS  NOT  SPECIFY  FI.OW  FIFLD  SURVEY  DATA  /  96. 

?  SX,  UHCARO  TITLFs.  14*)  97. 

P010  FORMAT  ( l HO .  SX.17M*  •  »  FRROR  •  »  *.  9P. 

1  SOHTAPE  TITLE  OOFS  NOT  SPECIFY  FLOW  FIELD  SURVFY  DATA  / 

?  SX.11HTAPF  T I Tl F=«  PSA*)  100. 

?0?0  FORMAT  ( 1 H 1 , 1 7H I NPUT  RUN  NUMHFRS  /  101. 

1  ?x . 1HI . IX, 7HRUN  NO.  /  ( 1 X . I P . F 1 0 . P ) )  10P. 

?010  FORMAT  ( 1 H l , 1 X  «  PSA* )  101. 

pn*0  FORMAT  ( 1  HO .  //  1X.7HRUN  NO. , 7X , pohparent  configuration.  /  in*. 

1  SX.I1.RX.7A*)  ins. 

POSO  FORM*T  ( 1  HO . 1 7X , 1*H***»*  FREE-STRFAM  CONDITIONS  ••••*,  /  in6. 

)  *X.*HMACH.*X .pHRT .6X ,?HTT , 7X. 1HP.7X , 1HT .7X . 1HV. 7X ,  107. 

P  1  HO , 7  X . SHRF /FT,  /  |0P. 

1  SX , ?HNO.*X , *HPSF A , 7X .6HRFG (R) . IX ,*HPSF A, IX .  109. 

*  6HOFG(R) ,?X ,6HFT/SFC.*x, 1HPSF.6X.6HFT (-1 )  /  110. 

S  PX.F6.3.6FR.] .F7.?.6H*10**6)  HI. 

?060  FORMAT  ( 1  HO , *  IX . ?7H»***«  FLOW  FIELD  DATA  •••••  /  11?. 

1  6X.?HXP.SX.PMYR.SX.?M/P.SX,3HALP.SX,1MSI6.  111. 

?  SX ,*HALPT ,SX . 1MRHI . 4X .*HVX/V. *X , *HVY/V.*X .  11*. 

1  4HV7/V.4X ,*HVI  /V.SX.PHML ,SX,*H0L/().*X,6HPTP/PT .  1JS. 

4  /  6X  .  PHIN.SX  .  ?MlN,  *X  .  PHIN.SX  .  1HDFG.SX  ,  IHOf  G.SX  .  116. 

5  3HHFG.6X  .  1HDF6  I  117. 

P070  FOPMAT  (  IX  ,F6.?.  ?  (  ?X  .F  S.P)  ,  1  (?X  ,FX  .?>  ,?X  ,f  7.?,  7  (  ?«  ,F6. 1)  )  1)«. 

POPO  FORMJT  ( inn. IX . 1 ?HA! PMA  ALPHA  /  SX.11HPAR  PROHF  /  ])Q. 

1  SX.inHOFG  OFG.  /  4X.FS.?«?X.FS.P)  IPO. 

r  IP1. 

r  HP. 

STOO  1?1. 

ENO  1 ?* . 
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SUROOUTINF  CONFIG  (O.  IfONF  .  ISTOwE  ) 

c 

c  OFCTPHFR  IHF  COM  IM.RaT  ION  C  OOF  Aflll  SF  T  UP  AN  AIPmAHFTIC 

C  ARRAY  FOB  IMF  CONFIGURATION  AMO  ImF  storf 

c 

c 

dimension  ir  (  1  a  i  .  icons  (  7)  .  IS  (  7 )  ,  [  s  1  uRf  (  u 

DATA  I  C  /4HF-4C .  4  H  ,  P  1  ,  ,  AH  T  •  (  S  «  4  Hp )  1 1 ,  .  4  H  P  )  F  .  , 

1  AM. Pf  .  AH  T  •  4H7||I  ,  4H(S7).  AMU  .  A  MF  . 

?  AH. POL.  AHF.PO.  AHP1F  .  4HTM.{.  4HS?)F.  AM. (SI. 

7  AH) F  / 

DATA  IS  /4HMK-R.  4H7.MF.  AH  3  •  HU .  4H7.AF.  AH7.AIJ. 

I  AHPROP .  AHF  / 

OATA  IBLANK  /AH  / 

C 

IF  ( R  .EO.  610.0)  St  Ql.fl 
HR  =  R 

MY=  MOD (HP. 10) 

HXs  HR/10 

C ASF  CONFIGURATION 
CASF  X=0 

IF  (HX  .NF.  0)  GO  TO  SO 
DO  AO  1=1,7 

I CONF ( I ) =  1  BLANK 
AO  CONTINUF 
GO  TO  ISO 

SO  GO  TO  (60. 70.R0.R0, 1 00. 1 10. 1?0. 130. 1  AO)  ,HX 

CAS F  X=1 

60  CONTINUF 

I  CONF  ( 1 )  =  ICO) 

no  ss  i -  p.7 

1  CONF  IDs  IBLANK 
6S  CONTINUF 
GO  TO  ISO 

CASF  X  =  ? 

70  CONTINUF 

I CONF ( 1  )  =  I C  I  |  ) 

I CONF (  P)  =  IC(S) 

DO  7S  1=  7,7 

I C  ONF  (  I  )  s  TPLANK 
7S  CONTINUF 
GO  TO  ISO 

CASF  X  -  7 

RO  CONT  INI  IF 

I CONF (11=  ir ( 1 ) 

TCONF(?)s  I r ( p  > 

I  CONF  (  7)  =  1C  (  n 
OO  RS  1=4.7 

I  C  ONF  (  I  )  =  Ini  AN« 
as  CONTINUE 
go  to  iso 
C 

C  CAST  X  s  4 

c 

RO  CONT INK) 

no  <is  1  =  1.7 


<>00. 
?oi . 

?Op. 

?03. 

?oa. 
pos. 
/Os. 
pn  i . 
?OR. 
?  OR. 
?1  0. 
PI  1  . 
PIP. 
?1  3. 

?  1  4  . 
PIS. 
P16. 
PI  7. 
PIR. 
?  1 B  • 

PPO. 
?Pl  . 

PPP. 
?P3. 
?P4. 
PPG. 
PP6. 
??7. 
PPR. 
PPR. 
P30. 
?71  . 
?7P. 
?73. 
?7a. 
?7S. 
P76. 
P77. 
P7R. 
?7R. 
P  4  0  , 
?A1. 
PA?. 
/A  7. 
P44  . 
PAS. 
?A6. 
PA  7. 

PAR. 

PAR. 

PS0. 
PS1 . 
?sp. 
?s  7 . 

PSA. 

?ss. 

?S6. 
PS  7  . 

P  SR. 
?SR. 

?sn . 
pm . 

P'-p. 
PS  ). 
p  s  4  , 
Pss. 
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ICONF  <  T  )  =  IC  (  I  )  ?SN. 

QS  CONTINUE  ?s7. 

TCONF  (4 )  =  IC  (  4)  ?*«. 

no  OS  1=  S.7  c*so. 

ICONF (11=  1  HI  ANK  ?70. 

OS  CONTINUE  ?f\» 

0,0  TO  ISO  ?7?. 

C  ?73. 

C  CASF  X=S  ?74. 

C  ?7S. 

100  CONTINUE  ??<>• 

no  ios  r=  i.4  ?77. 

ioonf<i>=  loin  ?7h. 

10S  CONTINUF  ?7<7. 

iroNFisi=  rc(oi  ?oo. 

rroNF(s)=  iron)  ?Bi. 

!OONF<7>=  I RL  ANK  ?«?. 

f.O  TO  ISO  ?R3. 

?R4. 

C  ASF  X=E> 

?RE>. 

110  CONTINUE  ?R7* 

DO  US  1=  1.3  ?««• 

ICONF  (  I  )  =IC  (  n  ?Rt>. 

11S  CONTINUF  ?*»<>. 

IC0NF(4)=  IC(S) 

ICONF(S)=  IC(O)  ?*»?• 

ICONF(N)=  I  C  < 1 1 )  ?Q3. 

I CONF (71-  T  RL  ANK  ?04 . 

00  TO  ISO 

?OS. 

CASF  X  =  7  ?“>7. 

?Ofl. 

1?0  CONTINUF  ?Qt,• 

IC0NF(1I=  Kill  300. 

IOONF<?)=  IC(1?)  301. 

I CONF  <31=  IC(?)  30?. 

ICONF (4) =  10(3)  30  3. 

I CONF ( S ) =  IC(S)  304. 

I CONF ( S )  =  1C(R)  30S. 

I CONF (71=  I C ( 1 3 )  3  OS . 

C,0  TO  ISO  307. 

30R. 

CASF  X=R  300. 

310. 

1 30  CONTINUE  31 1 . 

no  ]3S  1=  1.3  31?. 

ICONF  ( I ) =  I  C  (  I  )  313. 

13S  CONTINUE  314. 

ICONF ( 4  >  =  T C ( T  4  >  3 1 S . 

00  1 3S  I =  S, 7  31  N. 

ICONF (1)=  1  ML  AN*  317. 

1 3S  CONTINUF  3 1 M . 

SO  TO  ]S0  31R. 

C  3?0. 

C  CAST  «r<!  3?  1  . 

c  3??. 

140  CONTINUE  3? 3 . 

I CONF  (  I  I  =  1 C  (  I  I  3?4. 

ICONF  (?l  =  IC  (->1  3?s. 

ICONF ( I) =  ic ( IS!  3**. 
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trONF(4)=  iriiM  3?7. 

Iconf <  s  >  =  jr ( 1 7 1  . 

rroNF  (M=  ir*(  i  wi  i?6. 

ICONF  (  7)  =  l‘N  AN*  310. 

C  3  31. 

C  FNO  CASF  CONFIGURATION  33P. 

ISO  CONTINUE  333. 

334. 

335. 

CASE  STOPF  336. 

337. 

7STOPF(l>=  ISIl)  330. 

ISTOPF ( 3) =  IRLAN*  330. 

r,o  to  < 170. iso. ion,?no,?io.??o.?3o> ,  my  340. 

341  . 

CASF  Y=1  34?. 

343. 

170  CONTINUE  344. 

1ST  ORE I  ? )  =  IS(?I  34S. 

GO  TO  ?40  346. 

347. 

CASE  Y=?  34fl. 

340. 

180  CONTINUE  3S0 • 

1ST  ORE  (  ?  )  =  I  S  (  3  )  3S  1  . 

GO  TO  ?40  3S?. 

3S3. 

CASE  Y  =  3  3S4. 

3SS. 

IOO  CONTINUE  3S6. 

I  STORE  (?) =  I S ( 4  >  3S7 . 

GO  TO  ? 4 0  3S8. 

3SO. 

CASE  Y=4  360. 

36)  . 

?00  CONTINUE  36?. 

1ST  ORE ( ? ) =  ISIS)  363. 

GO  TO  ?4 0  364. 

36S  , 

C  ASF  Y  =  S  366. 

36  7  . 

?10  CONTINUE  368. 

ISTOPF (11=  T  S ( 6 )  360. 

1ST  ORE  I  ?  >  =  I S (  7 )  370. 

GO  TO  ?40  371. 

37?. 

C ASF  V*  6  377. 

??0  CONTINUE  374. 

1ST  OPE ( ? 1 =  ISI?)  3  7S , 

GO  TO  ?40  376. 

37  7. 

CASE  Yr7  178. 

?30  CONTINUF  370. 

(STORE  I?) =  I S ( 3 )  300. 

C  30  1  . 

C  FNO  CASE  STOPF  38? . 

c  343. 

?4 0  CONTINUF  384. 

C  IMS. 

RFTllPN  386 . 

FNO  387. 
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NOMENCLATURE 


The  following  list  contains  the  nomenclature  used  in 
the  report  except  for  that  used  in  the  model  component 
designation  and  in  the  data  tabulations.  That  nomenclature 
is  defined  in  lists  in  the  text  of  the  report. 


A  . 

z 


'R 


Cv  CLL 


C  ,  CLM 
m 


C  ,  CLN 
n 


C  ’N 

if*  1 


CY,  CY 


Engine  exhaust  exit  area,  7.4662x10  3 
square  feet,  model  scale 

Engine  inlet  area  used  in  data  reduction, 
1.705xl0-2  square  feet  model  scale 

Capture  ratio  calculated  using  Equation  (3) 

Axial-force  coefficient,  axial  force/q^S 

Store  rolling-moment  coefficient,  rolling 
moment/q^Sd 

Store  pitching-moment  coefficient,  pitching 
moment/q  Sd 

'  ^oo 

Store  yawing-moment  coefficient,  yawing 
moment/q^Sd 

Store  normal-force  coefficient,  normal 
force/q^S 

Store  side-force  coefficient,  side  force/q^S 

Store  diameter,  0.7  inch  model  scale 

Length  of  modified  afterbody  store,  5.723 
inches  model  scale 

Engine  exhaust  Mach  number 

Free-stream  Mach  number 

Engine  exhaust  mass  flow  rate  calculated  using 
Equation  (2) ,  lbm/sec 

Free-stream  total  pressure,  psfa 

Average  static  pressure  in  engine  exhaust,  psfa 
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NOMENCLATURE  (Contd. ) 


a 


i?t 


VY/V 


vz/v 


v 

:  •■>  -  p* 


"p 


Average  total  pressure  in  engine  exhaust,  psfa 
Free-stream  dynamic  pressure,  psf 
Free-stream  Reynolds  number  per  foot,  ft-1 
Radial  distance,  in 

Reference  area  which  is  store  frontal  area, 
0.385  square  inches  model  scale 

Free-stream  total  temperature,  °R 

Engine  exit  velocity  calculated  using 
Equation  (4),  ft/sec 

Free-stream  velocity,  ft/sec 

Local  Yp  velocity  component,  positive  in  the 
positive  Yp  direction,  divided  by  the  free- 
stream  velocity;  also  termed  sidewash  angle, 
rad 

Local  Zp  velocity  component,  positive  in  the 
negative  Zp  direction,  divided  by  the  free- 
stream  velocity;  also  termed  upwash  angle,  rad 

Coordinate  system  in  which  probe  static 
pressure  taps  are  located;  origin  shown  in 
Figure  5 

Coordinate  system  in  which  store  nose  is 
located;  origin  shown  in  Figure  5 

Assumed  flow  separation  location  on  store  body 
measured  from  nose,  in 

Store  vertical  position  relative  to  attached, 
or  carriage,  position;  equal  to  Zg/D 

Angle  of  attack  of  the  probe,  deg 

Angle  of  attack  of  the  store,  deg 

Store  rolling  angle,  positive  clockwise  looking 
upstream  and  measured  relative  to  fins  45 
degrees  from  the  vertical  and  horizontal,  deg 
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